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Zur Asymptotik der Zylinderfunktionen und ihrer Nullstellen 
| "Von WALTER FRANZ 


Die asymptotischen Darstellungen und der Verlauf der Nullstellenkurven der Zylinderfunktionen 
für beliebige komplexe Argumente und Indizes werden angegeben, ohne auf die schwierige Diskussion einzu- 
‚gehen, über welche Sattelpunkte jeweils die Integrationswege der Sommerfeldschen Integraldarstellung 

 verlaufen!). Statt dessen wird benützt, daß die Nullstellenkurven die Verzweigungsschnitte der asymptoti- 
schen Darstellungen sind; auf Grund dessen lassen sich die asymptotischen Formeln ebenso wie der Verlauf 
der Nullstellenkurven für beliebige Argumente und Indizes durch analytische Fortsetzung gewinnen. Für. 
nahezu reelle Indizes ist eine gesonderte Diskussion erforderlich, da für rationale bzw. ganze negative Indizes 
die Nullstellen der Hankelfunktionen in gewissen Zweigen nach dem Unendlichen, bzw. die der Bessel- 
funktionen in den Nullpunkt auswandern. 


Asymptotic representations and curves of zeros are given for Bessel functions of any complex arguments 
and indices. The difficult discussion as to which saddle points the path of integration in Sommerfelds 
integral representation has to pass through!) is avoided. Instead, the fact is used that the curves of the zeros 
are the branching sections of the asymptotic representations. Hence, asymptotic formulae may be obtained 
as well as the zero curves for any arguments and indices by analytic continuation. A special discussion 8 
required for indices that are almost real because in certain branches the zeros of the Hankel or Bessel 

- functions shift towards infinity or the origin if the index becomes rational or negative integer, respectively. 


Ina 6ecceJreBbIx PyHRINÜ HPOH3BOJIBHOTO KOMINIERCHOTO APTyMeHTa MU HUHAERKCA NAIOTCA 
ACHMTOTHYeCKHE IPeNCTaBAleHNnA U TOBeNEHHE KPUBbIX, COOTBETCBYIONIUX HYJAIAM ITUX DYHRIMH. 
3atpyiIHHTelbHbIe PAccysKNeHuA 0 TOM, yepe3 Kakne TuIepÖoAMYecKMe TOURH IPU 3TOM 
IPOXOAAT HYyTu MHTerPuPOBAHuA III HHTETPAIIBHOTO HPejcTaBaeHuA 3oMMepherbla, 3Tech 
usa0eramrtca!). BMecTo 3TOTO UCHOJIB3YeTCA TOT PART, UTO KPuBbIe Hyıleü ABAAIOTCA Pa3pesamHu 
Pa3BeTBJIeHuA ASIA ACHMIITOTUYCCKUX BbIpaskennfi. Ha ocHaBaHnn 3Toro ylaeTca HOTy4MTb 
ACHMUOTOTUYeCKME POPMyYYIbI U ONPEeNEeJIHTb MOBEeNEeHME KPMUBbIX, COOTBETCBYIOIHX HYJIAM, MI 
IPOW3BOJAbHBEIX APIyMeHTOB M WMHHNEKCOB IIPM TOMOIIM AHAJIMTUYECKOTO HPONOJBREHNUA. 
Cııyyali IIOyTU BeIlleCTBeHHbIX UHJEKCOB HY;RHAaeTCcA B 0CO00M paccy’KIeHUN, TAR KaR AA 
PallHOHAJIbHEIX, COOTBETCBEHHO IICJIEIX OTPHNATENIBHBIX UHJIEKCOB HyJIM PyHRUNÄ Xankesst Ha 
HEKOTOPbIX BETBAX YNANAIOTCH B ÖeCHOHEYHOCTB, COOTBETCTBEHHO HYJIU OeccelleBbIX PYHRIMH 
IIPNUÖJIMSKAIOTCH K Hayalıy KOOPAHHAT. 


1. Einleitung 


Die asymptotischen Entwicklungen der Zylinderfunktionen sowie die Lage ihrer Null- 
stellen kann man übersichtlich und mit verhältnismäßig geringem Aufwand an Mühe diskutieren, 
wenn. man beides gemeinsam betrachtet. Der Zusammenhang wird dadurch hergestellt, daß 
die Verzweigungsschnitte der asymptotischen Darstellungen identisch sind mit denjenigen 
Kurven, auf welchen die Nullstellen liegen (Nullstellenkurven). Auf diesen Kurven sind nämlich 


- asymptotisch die von zwei verschiedenen Sattelpunkten der SOoMMERPELDSchen Integraldar- 


stellung gelieferten Beiträge von gleichem Betrag, was einerseits dazu führt, daß beim Über- 
schreiten der Kurve der eine Sattelpunkt durch den anderen abgelöst wird, somit die asym- 
ptotische Darstellung wechselt, während andererseits beim Entlangwandern auf der Kurve mit 
der Veränderung der Phasendifferenz Punkte erreicht werden, auf welchen die beiden Beiträge 
sich kompensieren, bei welchen also Nullstellen liegen. Man kann für beliebige Argumente ‚die 
asymptotischen Darstellungen sowie die Nullstellenkurven der Zylinderfunktionen Z,(z) in der 
»-Ebene ermitteln, wenn man nur die asymptotische Darstellung für » > jz| kennt, und weiß, 
daß für positives z die Nullstellen der Bzsserfunktion auf der rellen Achse bei » < 2 liegen 
und für positiv imaginäres z die Nullstellen von H,’@) bei rein imaginären » mit | > |zl. 


1) G. N. Watson, Theory of Bessel Functions, Kapitel VIII. Cambridge 1952. F. Empe, ZAMM 19 


(1939), S. 101. P. BECKMANN u. W. Franz, ZAMM 37 (1957), S. 17. 
26 
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Die asymptotische Darstellung in der Umgebung der positiven z-Achse (verschwindendes 
Argument p von 2) ist 


7 a ? 
mw || —n Sin y— y Cos für > ee (1). 
24) os, ar 6m Cool E 
Dabei ist $ 
Cosy = - ee (2), 


y sowie das Wurzelvorzeichen werden eindeutig festgelegt durch 


P TT 
Re)>0; ImW+g<z; agGiny+el<z 


Daß J,(z) bzw. K,(z) für @ = 0 nur bei reellem bzw. rein imaginärem » verschwinden können, 
wobei »<z bzw. »| > |z| sein muf, folgt leicht aus der Besserschen Differentialgleichung 


in der Gestalt 


für Z=2J,, 3 

A 7 . (8). 

Re(y)<0; ImY)—4<5; arg Sin) tr —p<Z ( 
— H® \ 

für Z= _ H® 

. 


— 1 I) = (v? — 22) J (X) , zK(@d) =? +9)Kfd) :...Ü: 


Daraus ergibt sich zunächst 


LH) TE TO 


L 


und eine analoge Gleichung für K,. Wir setzen nun voraus, daß x reell und J,(z) = 0 ist. Weiter 
sei u zu » konjugiert komplex; dann muß auch J,(x) = 0 sein. Damit verschwindet die linke 
Seite von (5), also auch die rechte. Da der Integrand positiv ist, muß »® = u? sein. Es kommen 
also nur reelle oder rein imaginäre Indizes in Frage. Für J, scheiden rein imaginäre »-Werte aus, 
weil J,und J_, dann gleichzeitig verschwinden müßten, was ausgeschlossen ist, da für imaginäres » 
die WRonsKıI-Determinante von J, und J_, nicht verschwindet. Für K, kommen reelle »-Werte 
nicht in Frage, weil dafür der Koeffizient der rechten Seite von (4) positiv wird; da die Funktion 
im Unendlichen verschwindet, muß für ein x, welches größer ist als eine beliebige Schranke, 
Funktion und Ableitung entgegengesetztes Vorzeichen haben — dann kann aber wegen (4) die 
Ableitung und deshalb auch die Funktion das Vorzeichen nicht wechseln, wenn man zu kleineren 
Werten von x übergeht. Auch für rein imaginäre » kann das Vorzeichen erst wechseln, wenn 
der Koeffizient von (4) negativ, also |v| > |x| geworden ist. — J, verschwindet für positives » 
im Nullpunkt, also gibt es unterhalb jeder noch so kleinen Schranke ein positives x, für welches 
Funktion und Ableitung gleiches Vorzeichen haben, nach (4) kann bei größeren x die Ableitung 
und damit auch die Funktion ihr Vorzeichen erst wechseln, wenn der Koeffizient der rechten 
Seite negativ, also x > » geworden ist. 


2. Mögliche Nullstellenkurven und ihre Wanderung 


Wir beschränken uns im folgenden auf die einfachen Zylinderfunktionen 2 J,(z), — H®(z) 
und — H}”(z). Nullstellenkurven können asymptotisch dort auftreten, wo für zwei Zweige 
von Gl. (1) der Realteil des Exponenten gleich ist. Diese beiden Zweige können entweder dadurch 
auseinander hervorgehen, daß y an einem Punkt mxi gespiegelt wird (Verzweigungspunkte 
v»= + 2), oder durch Verschiebung von y um ein Vielfaches von 2rxi. Daß beide Beiträge vom 
gleichen Betrage sein sollen, liefert die Bedingungen 


Re [e? (Siny—- y—mnl)Cosyjj=D. Sn ee 
oder 


Re [i e’Cosy]=0; dh, Im) =: ne 
Man sieht, daß sich Kurven (6a) für verschiedenes m nur auf der reellen v-Achse schneiden 
können. — Wir setzen 

y—mnai=a-+iß; %,ß,rteell 


e 
E 
E z 
4 
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und zeichnen die Lösungskurven von (6a) bzw. (6b) in Bild 1 und Bild 2 auf; diese Gleichungen 
lauten, auf a, ö umgeschrieben: . 

Tan a — oa + ß Tan a tan ß 

Near nn (7a), 


TE Ahle Se er a re) 


tanpo= 


Jede einzelne Nullstelle ist asymptotisch durch eine feste Differenz zwischen den Exponen- 
ten zweier Zweige von (1) gekennzeichnet. Das Linienelement der Nullstellen-Wanderung erhält 
man durch Differentiation zu r 

Sin, d2)—= (1. se dp me een) 
bzw. N 
Ren DS 1 Br ee ee hi): 


Man sieht daraus, daß in Abhängigkeit von z die Nullstellen der ersten Art stetig in der 
v-Ebene wandern, die der zweiten Art asymptotisch ruhen. Ändert sich v, so wandern die Null- 
stellen der ersten Art nach (8a) in der z-Ebene stetig; um die zu gegebenem » gehörigen z-Werte 

der Art (8b) zu ermitteln, reicht die 
asymptotische Darstellung (1) — auch, 
wenn man die gesamte asymptotische 
Reihe hinzunimmt — nicht aus; hier- 
für muß man zusätzlich die Umlaufs- 
relationen der Zylinderfunktionen be- 
nützen. — Natürlich gelten die Glei- 
chungen (8) für die Bewegung der Null- 


Gleichung (6) die Nullstellenkurve wieder- 
gibt, bis also die asymptotische Dar- 
stellung (1) sich durch Hinzutreten eines 
weiteren Zweiges ändert, m. a. W. die 
betrachtete Nullstellenkurve von einer 
anderen geschnitten wird. 


Bild 1. Mögliche Nullstellenkurven nach Gl. (6a). tan p ist Parameter Bild 2. Reelle »-Achse in der y-Ebene. 
tan @ ist Parameter 


An einem Verzweigungspunkt » = + z entspringen drei Nullstellenkurven in den Richtun- 


gen VF =. Man erhält nämlich aus (6a) für y — mil < 1 die Bedingung Re [2 (y — mr ij] =0, 
d.h. — Y— mrüfz?>0. Daraus folgt.nun - 


ee Eee a a Re 
zZ 2 yF dar th-ma =; re bem AS TITz 


Damit ist die Behauptung bewiesen. Die drei bei + z entspringenden Nullstellenkurven der 
v-Ebene drehen sich also bei einer Veränderung von @ mit einem Drittel der Winkelgeschwindig- 
keit mit. 


3. Nullstellenkurven in der Index-Ebene 


Die Wanderung der beiden Ufer der Nullstellenkurven in der y-Ebene mit wanderndem z 
ist durch die Bedingungen (6a) und (6b) eindeutig vorgezeichnet. Die Lage der Kurven ist in- 
Bild 3 für eine Reihe positiver Werte von p angegeben und in die »-Ebene umgezeichnet. Für J, 
hat man von der Kurve bei p = 0 auszugehen, für H,” von @ = n/2, für HP” von g = —.n/2, 
um den Anschluß an den aus der Einleitung bekannten Fall zu finden. Wegen der Umlaufs- 
relationen der Bzsserfunktion wiederholt sich deren Nullstellenkurve, wenn @ um ein Viel- 
faches von x wächst. Die Nullstellenkurven der modifizierten Hankerfunktion K,(z) — d.h. 
von H®(i z), H®(—- iz) — enthalten das Mittelstück der reellen Achse, sobald |p| = 3/2 ist. 
Wächst |p| über diesen Wert hinaus, so nähern sich die Außenäste der Nullstellenkurve monoton 


der reellen Achse. 
26* 


stellen nur solange, wie die asymptotische 
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4. Asymptotische Darstellungen 


Außerhalb der Nullstellenkurven kann die asymptotische Darstellung (1), (2), (3) analytisch 
fortgesetzt werden. Für v» > |z| ist diese Fortsetzung der Darstellung von J, bereits durch 
Gl. (3) gegeben, ebenso für » Z 1,51 |z| die Darstellung der HAnkEer-Funktionen; denn diese 
Gebiete der positiven Achse bleiben von Nullstellen völlig frei. Von hier aus kann man nun in 
jedem Blatt der z-Ebene die Darstellung für sämtliche » gewinnen. Auf den Nullstellenkurven 
selbst gelten diese asymptotischen Ausdrücke nicht mehr, vielmehr hat man dort die Summe der 
beiden Zweige der Funktion (1) zu benützen, welche den beiden Ufern des Verzweigungsschnitts 
entsprechen. Die Darstellung durch zwei Zweige von (1) gilt jedoch nicht nur auf dem Verzwei- 
gungsschnitt selbst, sondern in einer sehr weiten Umgebung, in welcher dann jeweils einer der 
beiden Zweige gegen den anderen verschwindend klein wird. Diese Umgebung reicht bei den 
auf der reellen Achse gelegenen Verzweigungsschnitten bis ins positiv und negativ imaginäre 
Unendliche über und unter dem reellen Achsenstück (und sogar noch erheblich darüber hinaus). 


enkurven kann. die De Sen zwei - 
h n anderen i im Winkelabstand von 120° dort entspringenden 
rven. Es ist daher zur Aufstellung der asymptotischen Darstellungen 
die genaue Lage der Nullstellenkurven zu an wenn man einfach die 
ne ch zwei Summanden wählt. 


Gehen von einem Punkt der v-Ebene drei Nollstelleniureen aus, So aan dort drei 
B ige von (1); die asymptotische Darstellung ist die Summe aller dreier; sie kann auch außer- 
halb des Schnittpunktes verwendet werden in dem ganzen Gebiet, welches die Gültigkeitsbereiche 
‚der auf den drei einzelnen Nullstellenkurven gültigen Darstellungen gemeinsam haben. 


Etwas komplizierter sind die Verhältnisse auf dem Mittelstück der reellen Achse für die 
Hankeıfunktionen. Wird der Winkelwert 9 = 2 überschritten, dann treten — wie Abb. 3 
_ zeigt — zwei Zweige der asymptotischen Darstellung aus dem analytisch zugänglichen Bereich 

' aus; man würde deshalb eine falsche Darstellung erhalten, wenn man durch analytische Fort- 

setzung von außen her auf die Nullstellenreihe schließen würde. Man muß vielmehr diese ‚‚ver- 
 lorenen“ Beiträge bei wachsendem @ beibehalten. Beim Durchgang von 9 durch ganze Vielfache 
von rı treten dabei jeweils zwei weitere „‚verlorene“ Zweige hinzu, wie dies der Umlaufsrelation 


sin (m —1)vn Hol) — er sinmvn Ho) 


sinv sc sinvn. 


EN emzi 2) a3 


entspricht; es gilt nämlich 


sin mv st 


— em 1)iva 4 ei(im—3)ivn ie an a eim—1)ra | 
sın v 7 


Für nahezu reelle Werte von z sind dabei Sattelpunkte zu berücksichtigen, welche z.T. 
auf Grund von (6a), z. T. auf Grund von (6b) gleiche Höhe haben. Daher wandern die Null- 
stellen bei einer Veränderung von z, wenn auch in einer von (8) verschiedenen Weise. 

Nicht erfaßt wird von den bisherigen Formeln derjenige Teil der v»-Ebene, welcher von den E 
_ Punkten y= +z einen Abstand von der Größenordnung z!/? oder kleiner besitzt. Hierfür gilt 
die Darstellung (1) nicht. Man muß statt dessen die Formel von NIcHoLson verwenden, welche 
von ScHöBE?) verallgemeinert und verbessert wurde, und geschrieben werden kann?) 


E Z „|wo-2- 3 : 2 0— a 11/8@ =) 
E AR "7 — = z RAD: 
ur. ld, Fk ee 


0 bedeutet den Differentialoperator 0/0», und die Phase des Arguments der Funktion Kıjs ist 
festzulegen durch die Vorschrift 


(v — 2)? 5 : tt 
arg Ba re 0 gegenüber Verzweigungsschnitt. 


Führt man die Differentiationen aus und entwickelt bis zur Größenordnung 283 (wobei 
man dem Argument von Kıj; keine Größenordnung zuschreibt), so ergibt sich 


oe a (PR + 2 — ) N RT) 


mit den Abkürzungen 
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2) W. SCHöBE, Ann. Math. 1 (1948), p. 220—232. 
> 3) W. Franz und R. GaLLe, Z. Naturforsch. 10a (1955), S. 374. 
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5. Nullstellenkurven in der z-Ebene 


Von den beiden Nullstellen-Bedingungen (6) ist nur (6a) genau genug, um die Nullstellen- 
reihe in der z-Ebene anzugeben. Da wir jetzt v festhalten, führen wir dessen Phase y in (6a) ein, 


und erhalten j 
Re [e’ (Tany—  — mazi))] = 0 
Sin2#«—a(Cos2« +coS2f)t -- «re... . (13). 


ANY sn 28 —B(Cos2a + cos2) 


dan: 


Den Verlauf dieser Kurven in der «-f-Ebene zeigt Bild 4. Die in der Umgebung von 
+3rij/2; +5mi/2 von den Kurven tan y = 0 eingeschlossenen endlichen Gebiete der a-ß- 
Ebene lassen wir aus, weil wir aus der 
Diskussion von Abschn. 2 und 3 ent- 
nehmen können, daß sie für Null- 
stellenkurven vom Typus (6a) nicht 
in Frage kommen. In diese Gebiete 
fallen vielmehr die auf der reellen 


gelegenen Nullstellen, deren Lage 
asymptotisch von z nicht abhängt. 


Es ist nun leicht, auf Grund von 
Abschn. 3 die für die einzelnen Zy- 
linderfunktionen zutreffenden Null- 
stellenreihen der Art(6a) auszuwählen. 
Das Ergebnis zeigt für einige Lagen 
von » Bild5 für die BzsseLfunktio- 
nen, Bild 6 für die HaAnkerfunktionen. 
Dabei müssen wir zunächst annehmen, 
daß die reellen »-Werte sich von den 
Nullstellen der Art (6b) fernhalten. 


In der y-Ebene sind die un- 
geraden Vielfachen von izz/2 nunmehr 
Verzweigungspunkte (z = 0), sie wer- 
den daher stets von einem Verzwei- 
gungsschnitt durchlaufen. Für die 
Bessefunktion zerlegt dieser Schnitt 
die gesamte y-Ebene in zwei Hälften, 
von denen nur die rechte verwendet 
wird. Bei H}’(z) kommt ebenfalls nur 
eine halbe y-Ebene zur Verwendung; 
diese verzweigt sich bei ir/2 in zwei 
Rıemanssche Blätter, wobei die asym- 
ptotische Darstellung (1) in jedem Blatt 

Bild 4. Mögliche Nullstellenkurven nach GI. (13). tan y ist Parameter bis zu einer bei in/2 entspringenden 

Kr N Linie gilt, welche sich periodisch um 
die imaginäre Achse schlängelt. Die beiden Blätter der Rrımmansschen Fläche überdecken sich 
dabei teilweise in lappenförmigen Ausbuchtungen. In der Figur sind die überdeckten Lappen der 
linken Fläche punktiert gezeichnet. Die Phase von z, welche den einzelnen Punkten der y-Ebene 
zuzuordnen ist, wird dadurch festgelegt, daß der Punkt y = 0 bzw. rn der Halbebenen-Begrenzung 
immer im Gebiet O0 <9< x liegt. ' 5 


v-Achse im Mittelstück (bei festem 2) 


ce u ee 
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| 


W. Franz, Zur Asymptotik der Zylinderfunktionen und ihrer Nullstellen 393 


2 @) en 
v2) 
=»p +P 
7 ee un ® 
u 
a) b) 


e) 


Bild 5a—e. Nullstellenkurven der BESSELfunktion in der Argumentebene für verschiedene Phasen y des Index » 


6. Besselfunktionen von nahezu negativem Index 


Nähert sich » der negativen reellen Achse an, so bleiben die in Bild 5 eingezeichneten Äste, 
welche sich dem Außenteil der reellen z-Achse nähern, richtig; denn sie genügen der Bedingung 
(6a). Die bogenförmigen Mittelstücke dagegen werden zweifelhaft, weil sie denjenigen Punkten 
von Bild 3 entsprechen, in welchen die beiden Nullstellenkurven b aus + z, —z und oo endigen, 
wo also die Bedingung (6a) durch (6b) abgelöst wird. Um dieses Übergangsgebiet zu erfassen, 
muß man die drei asymptotischen Darstellungen der angrenzenden Gebiete addieren. Es genügt 
dabei, den Fall In) >0 und O<e<n zu betrachten, da mit J,(z) = 0 auch J,(— 2) = 0 
und J,«(z*) =0. Nun gilt rechts von dem Schnittpunkt auf der gesamten reellen Achse und 
in ihrer Umgebung, wie man durch analytische Fortsetzung von (1) erkennt; 


a N Sure I m a AA). 


Wir haben also links vom Schnittpunkt zu der dort durch Fortsetzung von (1) erhaltenen, 
aus zwei Summanden bestehenden Darstellung noch die asymptotische Darstellung von e'”” J_, 
zu addieren. Man sieht aber leicht, daß man durch diese Hinzufügung auch die in Gl. (6b) fehlende 
Abhängigkeit der auf der negativen v-Achse gelegenen Nullstellen von z richtig erfaßt. Nach (6b) 
sollten die Nullstellen nämlich bei ganzzahligen »-Werten liegen. Für ganzzahlige Werte von » 
gilt aber Gl. (14) streng. Fügt man also den (sonst auf der negativen v-Achse asymptotisch zu 
vernachlässigenden) Anteil (14) der asymptotischen Darstellung hinzu, so erhält man den rich- 


E \ 170 e ea E 4 
bik der 2 ylinderfunktionen und ihrer Nullstellen“ ie: 


eh 


Nullstelle in Abhängigkeit von z berechnen. Dafür ergibt sich dann die Gleichung 


_ Die linke Seite ist asymptotisch für Im ) > 0 vom Betrage 1. Daß auch die rechte Seite 
3etrag 1 habe, ist die Bedingung, welche zu dem Zweieck von Bild 5e führte. Daher liegt 
‚diejenigen v»-Werte, für welche die linke Seite von (15) den Beträg 1 hat, die NStKurve auf 
dem Zweieck; ist sie von kleinerem Betrage, so liegt die NStKurve innerhalb des Zweiecks, ist 
sie dagegen merklich größer als 1, dann gibt es 
keine NSt in der Umgebung des Zweiecks. In 
_ der negativen »-Achse hat man somit nur dann 
NSt auf bzw. im Inneren des Zweiecks der z-Ebene 
_ zu erwarten, wenn die linke Seite von (15) dem Be- 
_ trage nach gleich oder kleiner als 1 ist; dadurch 
_ wird eine gewisse Umgebung der ganzzahligen 
 v-Werte gefordert, welche in Bild 7 schraffiert 
eingezeichnet ist. N 


- Denver aan, für Do (15). 


Ag 


- Die Lage der Nullstellen innerhalb des Zwei- 
_ ecks erkennt man am besten, wenn man Gl. (15) £ 
‚durch z und » ausdrückt: 
IE —2» 
er 2sinvyam( —— — | .eP=? (16). 
Er: 2, 


Die Wurzel ist mit positivem Realteil zu 
nehmen. Für » nahezu —.n ergibt dies 


R ei an { 
y _ n et De A eln— 2° (17). 


an Ym—z2+n 


7 


Bild 7. Gebiete der v-Ebene, in welchen Nullstellen gegen den 
Nullpunkt der z-Ebene wandern 


 zmuß daher in der Nähe des Nullpunktes liegen, 
und man hat genau 2n Lösungen, welche sich 
um die 2n-ten Einheitswurzeln unterscheiden, 

also ‚auf nem Kreis um den Nullpunkt gleich- Bild 8. Wanderung der Nullstellenkurven der BESSELfunktion 
mäßig verteilt sind. Nach der genaueren For- bei Annäherung des Index an eine negative ganze Zahl 
mel (16) hat man ebenfalls 2 n Nullstellen, welche 

jedoch nicht gleichmäßig auf einem Kreise angeordnet sind. Nähert sich » aus der positiv imagi- 
nären Halbebene dem Wert —n an, so verlassen 2n Nullstellen das Zweieck von Bild 5e und 
wandern in den Nullpunkt, wenn » bei —n einläuft. Dort kompensieren sie die Unendlichkeit 
z’ von J_, und schaffen stattdessen die n-fache Nullstelle von J_,„ = Jun. Dieser Vorgang ist 
in Bild 8 dargestellt. Lage und Anzahl der Nullstellen entspricht dem Hurwırzschen Theorem. 


7. Hankelfunktion von nahezu reellem Index - 


Nähert sich der Index der HAnkerfunktion der reellen Achse, dann bleibt die Diskussion 
von Abschnitt 4 nur für die zum Nullpunkt verlaufende Nullstellenkurve gültig, und auch für 
sie nur in ihren ersten beiden Quadranten (von z= + v ab gerechnet); denn nur in diesem Fall 
rührt die Nullstellenbedingung von (6a) her. In allen anderen Fällen befindet man sich in dem 
Gebiet, in welchem entsprechend Gl. (9) eine größere Anzahl von Zweigen der Funktion (1) 
zur asymptotischen Darstellung beiträgt. Wir benützen deshalb die Umlaufsrelation (9), um 
- die HankeLfunktion Hf”(z) auf ihre Darstellung in der Halbebene 0 <g <r zurückzuführen. 
Wir wollen dieses Gebiet das nullte Blatt der Rremansschen Fläche nennen und allgemein das 
Gebiet mn <p < (m— 1) das m-te Blatt. Wir setzen für einen z-Wert des m-ten Blattes: 


ER 


z, gehört dann dem nullten Blatt an. Für »-Werte, welche von der positiven Achse nur weni® 


entfernt sind (Abstand < v!/?), gilt — wie man aus Bild 3 entnehmen kann — im gesamten nulltel 
Blatt der z-Ebene die Darstellung (1), welche man (zugeschnitten auf das Gebiet |z,] > ») in der 
Form schreiben kann 


Ten Ai 
IN Zu) m Vs . 
( 0) De Vz} Fr 


D.Ssargı, Sr; 0. <arg (Ya —»v) <n 


ee I he BE Kl a A ae 


V3— +iv 


iVz2 v2-+-vlog 
e w. 


Se ee, 


n asym ptotischen Wert von J, bei ganzzahligem », und kann die in der Nähe dieses Punktes 


Fe = # a Te rn 


396 W. Franz, Zur Asymptotik der Zylinderfunktionen und ihrer Nullstellen 


ie zwei ion gilt für j iti i ten Blatt die Darstellung, 

Für die zweite Hankerfunktion gilt für jedes positive » im ers 
welche auf der Nullstellenreihe h, zu nehmen ist, nämlich die Summe der beiden Zweige, a: 
man von Gl. (1) aus durch Umlauf beiderseits um den Punkt z = v erhält. In dieser Weise folgt 


NEE TEEN a Y 2 —v? +ir Ne 

B“ Ho ya \ 2 a gr + ! 0) 
; m — Ay NZ, =, A A . 
Bl k ya — 


0<argau<r; 0<ar(Y2 —) <a 


i i i i i klein gegen 1 
Liegt z, so nahe an der reellen Achse, daß die Exponentialfunktion nicht mehr kle‘ Y 
ist, so Ei der Term — H“ auf der rechten Seite zu streichen. — Hieraus folgt als Bedingung für 
die Lage der Nullstellen 


Bu 2y 

1 2 2 j mi :Y. . 

3 ee 1v 3iVz2 —»: 

. eimvz u Vz 3 r ‚E13 F ar Te a a 21) 
sinv zz 29 

für komplexe z und 


inmen |_1 0.00. 
‚sin (m— I)v 


für reelle z. Da die rechte Seite von (21) in der positiv imaginären Halbebene kleiner als 1 ist, 
findet man die Nullstellenreihe nur oberhalb der reellen z,-Achse, wenn auch die linke Seite 
kleiner als 1 bleibt. Der asymptotische Wert der linken Seite für großen Imaginärteil von v 
nämlich i e@m—Drz, ]iefert genau die in Abschnitt 4 diskutierte Nullstellenkurve, welche für 
zwei »-Werte knapp oberhalb und unterhalb der reellen Achse in Bild 9 angegeben ist. Für reelle » 


+v 


b) 


Bild 9. Veränderung der Nullstellenkurven von Z/(l%z), wenn » die reelle Achse passiert 


schließen sich die beiden Mittelteile der Kurven von Bild 9a zu dem Zweieck von Bild 5e zusammen, 
um sich sodann von den Punkten + » zu lösen und eine Spirale zu bilden. Die rechte Seite von (21) 
nimmt in einem größeren Gebiete stark ab, wenn man sich von der Nullstellenreihe nach oben 
bzw. rechts entfernt (Rechtsschwenkung der Nullstellenkurve um den Punkt »), und sie nimmt 
zu, wenn man sich dem Zweieck bzw. der negativen Achse annähert. Liegt » in der Nähe einer 
Nullstelle der linken Seite (d.h. » unganz, m» ganz), so wandern die gesamten Nullstellen nach 
dem positiv imaginär Unendlichen der z-Ebene aus; dies hängt damit zusammen, daß für den 
genannten rationalen Wert von » die HAnKkeLfunktion im m-ten Blatt bis auf das Vorzeichen 
der im nullten Blatt gleicht, und daher von Nullstellen frei ist. — Nähert man sich einem reellen 
Wert von », für welche die linke Seite von (21) größer als 1 ist (größer als der bei ganzzahligem » 
auftretende Wert vom Betrage m kann sie nicht werden), so wird beim Erreichen des Betrages 
entsprechend (21a) das Zweieck bzw. die reelle Achse erreicht und sodann ein geringes Stück 
unterschritten. Nach dem Verlassen der reellen »-Werte in Richtung auf negativen Imaginärteil 
wird aber auf jeden Fall die linke Seite von (21) in der Umgebung des linken Außenabschnitts 
der reellen Achse wieder klein, und dieser Teil der Nullstellenreihe wandert nach oben, ent- 
sprechend Bild 9b. 


8. Sehlußbemerkung 


Die angegebene Untersuchung der Nullstellenreihen an Hand der asymptotischen Dar- 
stellungen liefert eine schöne Übersicht über die Nullstellen der Zylinderfunktionen, und gibt 
im übrigen Hinweise auf einige Sätze, deren exakte Gültigkeitsgrenzen für nicht sehr große 
Werte von z und » untersucht werden sollten. So wäre es z. B. interessant zu wissen, von welchem 
Werte »/|z| an Hankerfunktionen von positivem Index keine Nullstellen besitzen können; 


asymptotisch ergibt sich hierfür 1,51. 


Manuskripteingang: 6. 6. 1959 


Anschrift: Prof. Dr. W. Franz, Hamburg 36, Jungiusstr. 11 
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E Be m z ZAMM 40 (1960) Em 9, Seite BE 
B On the Solution of an Integral Equation Appearing 

a In the Delta-Wing Theory 

| By MILoMIR M. STanısıc!) | 


 Dedicated to Prof. Ludwig Prandlf for his enormous help during my studies 


‚Die Integralgleichung für den nichtstationären Fluß wm einen Delta-Flügel bei Überschallgeschwindig- 
keit wird untersucht und mittels eines neuartigen Approximationsverfahrens gelöst, das hier entwickelt wird. 


Er -; { Die Lösung der Integralgleichung repräsentiert die Quellenfunktion für den nichtstationären perfekten drei- 


dimensionalen Fluß. Die vorliegende Arbeit gibt die Lösung für den Fall der Überschallgeschwindigkeit, 
doch ergeben sich keine Schwierigkeiten für die Lösung außerhalb dieses Gebietes. Das neue Näherungsver- 
fahren basiert auf der Störung des Integrationsgebietes. Dadurch wird die gegebene Integralgleichung redu- 
ziert auf die bekannte Gleichung der Trragflächentheorie von Fuchs, Hopf und Seewald. - 


; Allgemein gesprochen bedeutet wohl die Entwicklung dieser T'heorie eine neue sich lohnende Methode auf 
diesem Gebiet, da sie einen Beitrag liefert zur theoretischen Erfassung der Hauptjaktoren, die beim ‚‚flutter‘- 
P Phänomen der Raketen eine Rolle spielen. 


The integral equation appearing in the formulation of the non-sieady motion of the flow around a 
delta-wing in supersonic flight has been studied and solved using a new approximate technique developed 
in this paper. The solution of this integral equation represents the source function for the three-dimensional 
non-steady perfect flow. In this paper the solution for transonic case is obtained, but no difficulties arise 
for the solution in the domain outside this region. The new technique developed in this paper for the solution 
of this integral equation is based on the perturbation of the domain of integration. This method reduces the 
rn ur under consideration to known equation in the aerofoil theory developed by Fuchs, Hopf 
and Seewald. 


Generally speaking, the development of this theory is believed to be one of the worthwhile new techniques 
in this field, which is a contribution to the major factors arising in the theory of the flutter phenomenon in 
rocket motion. ; 


HHrterpasaIpHoe YPaBHEeHNe, BCTpeyalmımeecH B VOPMYJIHPOBKE HEYCTAHOBUBIUETOCH ABU- 
SREHUA TeYEeHHA BOKPYT ACJIBTA00PAa3HOTO KpblIA IPH TOJAETEe CO CBEPX3BYKOBOH CKOPOCTBIO 
ÖbIJIO U3y4YECHO H pellieHo IIPH IIOMOINM HOBOTO IPHÖJNHFKEHHOTO METOAA, H3IIOFKEHHOTO B HACTON- 
ımeii pa6Oote. Peimenue 9TOrO UHTETPAABHOTO YPABHEHNA IIPencTaBlneT Co0oi PYHRUHIO NC- 
TOYHYURA MIA TPEXMEePHOTO HeyCTaHoBHBLIeTocH Teyennun. B HacTonmmeii pa6oTe BBIBOAHTCH 
peIueHne MIA CBEPX3BYKOBOTO CAIYYaA, HO U MIA pelieHuA BHe 3TOU O0JJACTU HU KARUX 34- 
TPyAHeHRü He 1oA1y4yaerca. HoBbliü MeToA, usJIaraeMmblÜ B HACTOAIMeH pabore In peimenunst 
HUHTETPaJIbHOTO YPABHEeHNSI OCHOBAH HA NPEOÖPA30BAHHN OÖJTACTU HHTETPHUPOBAHUA. HTOT 
IIPHeM CBOAUT PAaccMaTpuBaeMoe HHTETPAAIbBHOe YpaBHeHHEe K YPABHEHHIO, M3BEeCTHOMY N3 
TeOPHH BO31yIUHOU PosIbLu, pasBuToü DyKcoM,Toup$boM HU BeBalıbaoMmM. 


Booöme TOBopA, PasBuTue 9Toli TeOPHUH CUUTAeTCH ONHUM N3 CAMBIX MOYHBIX METOAOB B 
9TOH OÖSIACTH, YUTO CIIOCOÖCTBYET N3YUEHNMIO HANOOSICE BAHKHLIX BOIPOCOB, BOSHUKAIOIIUX B 
Teopun ÖHeHnü IIPu HBESKCHHH PareT. 


Introduetion 


Flutter is a fundamental consideration in the design of a transonic or supersonic aircraft. 
The current obstacle to the performance of an accurate mathematical analysis of the flutter 
problem in the transonic and supersonic range is the lack of an adequate theory of the aero- 
dynamic forces on an oscillating wing. The current applications of the delta wing in aircraft 
and anti-aircraft rockets are likely to continue. Thus, the distribution of aerodynamic forces 
on the delta wing is a basic problem. Any progress in the solution of this problem should be 
relevant to the study of re-entry vehicles. i 

In [1] Köüssner formulated the theory of flutter in 1929. Results obtained before 1950 
were limited to solutions of two dimensional non-steady state problems. The solution for the 
source function for non-steady state flow in three dimensions has not appeared although this 
field has not been neglected in recent years [2], [3], [4]. The difficulty lies in the solution ofan 
integral equation which results from the boundary conditions on the lifting surface [5]. 

The author of this paper is very thankful to the editor of this journal for recommending 
the new book [4] by J. W. Mızes, which appeared after this paper was submitted for publication. 
The theory of non-steady motion of the wing in supersonic flow has ‚been developed inde- 
pendently by the author [6] in the time from 1950—1958 during his association with the Mathe- 
matical Department of the Illinois Institute of Technology. J. W. MıLEs obtained the potential 
function by the use of the FOURIER transforms [7]. Therefore, all the integrals used in formulation 


1) Associate Professor of Engineering Sciences, Purdue University, Lafayette, Ind. 
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of the Integral Equation contain higher order singularities, which are inconvenient to evaluate 
for the aerodynamic force. J. W. Mıtes [4] did not solve integral equation resulting from the 
boundary conditions; he has suggested the method of Warkıss, C. E. and Berman, J. [8], [9]. 


An inspection of these solutions shows that the techniques for the steady and non-steady 
state solutions are the same, i. e., the region of integration is divided into strips such that the 
force can be assumed to act on each strip; once this is done, the superposition principle is 
applied to obtain the complete solution by the summation of the solutions for each strip. This 
result is due to HEASELT, Lomax, and Jones [2]. Since this approach cannot be used for non-steady 
flow, it is best regarded as an engineering approach. Thus, [8], [9], overcome the difficulties 
caused by the domain, which is the main obstacle in the solution of this integral equation. Here, 
a contribution is made in this direction. ‚The strip theory“ is rejected as a theory for unsteady 
flow, in the search for a direct solution of the integral equation with respect to the domain under 
consideration. A new technique in the theory of integral equations is required for the solution 
of the integral equation appearing in the theory of non-steady motion of the wing. 

An exact solution of this integral equation is not known. Communication by the author 
of this paper with $. CHANDRASEAKAR, L. COLLATZ, W. SCHMEIDLER, F. HILDEBRAND, S. S. SHU 
and M. GoLoMmB, has confirmed the opinion that the exact solution should not be attempted and 
an approximate solution should be obtained. A transform in form of hyperbolic functions can 
reduce the equation under consideration to a new form with unit kernel. However, the region 
of integration becomes more complicated. In other words, when we relax the kernel, we compli- 
cate the region of integration; the net result is that the difficulties increase. In that which follows, 
an approximate solution will be developed. 


Formulation of the problem 


Using linearized theory, the motion of the compressible non-steady flow is governed by 
equation 


1 
UM) Dat Dt Dr zB UBa+B)-=0 .......(), 


ö rk : 19: 
where ® is the velocity potential, M = 7 the MacH number, U is the velocity of the wing, 


and a is velocity of sound. The solution of this equation for a delta wing in supersonic flow 
can be obtained using the VOLTERRA-GREEN’S theorem [6] in form of acceleration potential; 
namely, 
X Bea e-ieot rn) 2 ’ 
> Fe ee | Yyı, 0) O&, y, 2; 2, y1, O0) dx, dy, - - -» . . (2). 
II* " 
where A(t,, 1, 0) is the source function for three-dimensional non-steady flow and 
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1 Iere A is the eigenvalue of the problem and ß® = M®— 1. If the position of the wing = 
any time is given by e2 


. 8). 


zuNuWeti 2.2... 2) 


where 9x, y) is a prescribed function regular everywhere in the region under consideration 
and © is the frequency of the vibrating wing, then, according to the linearized theory the con- 
ditions that the normal component of the velocity of the wing must be equal to the normal 
component of the velocity of the flow lead to the equation: Pr 


w(z, y, 2, D,o = ei [U gs, y)—iw Kr. ee 
>. > Y rm > > ‘ rare ır’# ” “ 
Using PrAnprL method of acceleration potential, we have that 
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Equations (2), 8), OÖ): and (6) after obvious substitutions in which care is taken in treatment 
of singularities during differentiation (and, if the resulting functions under consideration are 
expanded as a power series in the frequency) lead to 


2 *“ N. N 
| 
| | D) Aula, u) “)| 3) Ka, y; 2,3) ar) da, dyı. . (M. 
n= n= n=0 
- ]II* . 


The asterisk before integral sign shows that integration has to be considered in the HADAMARD 
sense. j 


Fig.1. Coordinate system of the wing 


From Equation (7) the system of integral equations result: 


; 1 ; Aylzı, Yı) @ — 2) dx, dy, 
Go(%, y) el ee - EEE (5) 
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The integral equation (8) represents the solution for steady flow given in [2]. This solution 
for steady case is obtained using ‚strip theory“. However, the solution of equation (9) with 
respect to A,(z,, y) represents the coefficient of A in the expansion of the source function. 
All other integral equations resulting from equation (7) are of the same type as equation (9), 
only they have other known functions on the left side. Therefore, if equation (9) can be solved, 
the complete system of integral equations can be solved. 


Henceforth, the attention is focused on the solution of the integral equation (9). The 
domain of integration is bounded by the leading edge of the wing and intersection of the for- 
cone with the wing surface (Fig. 2 and 3.) 


It must be pointed out that equation (9) is a singular integral equation with an unsymmetric 
kernel. The main difficulties in the solution of this integral equation appear in the complicated 
domain of integration and also in the complexity of the integration in the HADAMARD sense. 
However, the boundaries of the region of integration are linear functions, which suggests that 
a solution can be found. Hence, a promising fact for its solution exists. Since an exact solution 
of this integral equation is hard to find, an approximate solution is found. It can be mentioned 
that atransformin theform of hyperbolie functions can be used which reduces the integral equation 
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to one with a kernel which is equal to 1. Although the kern 


el becomes trivial, integration becomes 


more complicated and the diffieulties are increased. The new method used in this paper, con- 


sisting of the perturbation of the domain of integration, 
equation in the airfoil theory developed by Fuchs, Horr, and Sezwaro [10]. 


reduces this equation to the known 


Fig.2. General view of the geometry of wing plane in eones A and 7’ 


Solution of the integral equation 
The integral equation (9) can be written as 


5 1 [f Alıı yı) — yPT'2 
ee ee | 
s i ji ee ß ey dx, dy, ER.) 
Fir 
After expansion, equation (10) becomes 
H&y) = — - | Ayla Yı) IE % r Er +..-|dr,d 
je GT Te us 
Note that 
Pu—-y)<a—ı. 
Therefore, the series under integral sign converges. 
However, 
Gay) = ray). (12) 
, N 
Alu Yı) = 3 Pt? Ara Yı) (13) 
and the domain of integralion can be written as: 
N 
I*= 3 ß"D, (14) 


2 R =() 
Thus equation (11) becomes: } 


| Rn N 
ee y) = — | | Alla, | | oil 
= & 1) (y—yı)" 5 a) +++. | dx, dy, 


FL Dn 
Nn—=U 


. (15). 


#2 
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Therefore, from equation (15) it follows: 


NE.) 
TE ER Aa Alt Yı) 
1(z, y) = J pt EEE EU (16), 
z p®D, 
n=0 
At ‘ - 
= er) gay, = BEZ NEE a (17). 
@—1) Yy— yı)? 
N 
E B"D, Se) 
n=0 n=0 = 


ee En En Ta ie) a re er ee Va 


wu EA ar uni ae 


Iix,.y) 
WING PLANE 


- 
X 


Fig.3. Geometry ofregion ofintegration, IZ*, of the wing 


Evidently 
A,lt, Yı) = P? Asa, Yı) + RAR. Y) +R Aa Yı)tr. 0. (18). 
Hence, the system of integral equations (16), (17)..., can be solved if the solution of the 
integral equation (16) can be found. 
It must be pointed out that equation (16) has to be considered over the entire domain 
of integration for all values of , namely, equation (16) can be written as 


1 Ad(ltı yı) 1 Aa, Yı) 1 [ Ast, Yı) 
Gil, = — — —— ——- dr, d => did ——— ne dx, du, --- 
1(&, 9) = r Yı — ı @Yı a) y—yı)? - Yı 
Wr 


y— y—yı) 
D, ; PD, 
(19). 
In order to solve equation 2% it is necessary to solve the equation: 
Ada yı) q 
GI(x, . d er 8 NO BEER UNSERE 20 
@y) = ee 2, dy, (20) 


for A$(z,, yı)-. Then using iteration process, Aö(x,, y,) can be found to satisfy equation (19). 
In this paper, the solution for the transonic region will be considered only, since the solution 
for the entire supersonic region can be evaluated by taking a suffieient number of terms of (19). 
For the transonic region M approaches unity through values greater than one; i. e., for the 


transonic region 
ne. 2 raten 2h). 
B>0 
27 


ga-a+s er 


have been rotated about the point Ckx, en in such a way ae 


wi jr 


it follows that the boundaries of integration are: 
y=m2» y=-—ma 


Fig.4. Region ofintegration in the case of transonic regime 


Equation (22) is a new type of integral equation whose theory does not appear in the litera- 
ture. It is a singular integral equation with a variable domain of integration. Both linear and 
‚non-linear transformations are used to solve this equation. Theregion ofintegration (see in Figure 4) 
by triangle OAB can be mapped into region O’A’B’ (see Figure 5) by the linear transformations 


y=hS+tÄN 
ie anf Fer an 5), 
yyrdrtın \i 
where a,, 1, Y1, 9, are coefficients which are determined from the conditions: y 
OA=0BRz1l 1... 000 Pe (26). 
Furthermore, 
0>0 a 
A— 4’ 
B>B 


where here and subsequently the arrows mean „maps into“. From equation (25), it follows 
that: 


4% — Pıy 

a) & a EIER 
a 5, —Yı ßı 

hm AT 
ER a Ih —Nbı Pı 


a a u a nen denne rn 


Diit.ue 7% RT sie i 


reden 


Fig.5. en Hl inlegreklen Kom Ei. 4 
g- PA Yıtmaz 
s 2mz ’ 1 gmz 


” 


- 


we 


ons (31) and to , 


ET a Te 


00, > 0% 0) 


Au,nmz) — Aıı,o 
N De 1 EHER (34) 
m equation (27a) it follows that 
Er 83). 
e that nt Ag,—nz) satisfies equation (35), namely 
a ee (36). 


Ba,nn> Bo, 1) 


| Re Tr, 6 RA re BT - 
From equations (27b)- and 137), it follows that the point Ba,m.) salisfies: 
Bern Km—ylz=zahdı—Yıhı- --»--- ern. rnn (38). 
Equations (36) and (38) lead to 
| Br F z Be a 0 en (39). 
- However, equations (30), (33), and (39) lead to 
4 f a en (40). 
Be igeton 60) er and (40), in nam with (36) give: 
\ ee ” (41). 


The region of integration as shown in a can be ner into 
by using the mappings | - 


= 
je 


E e 
—m | ee Br 


Fig.6. Transformed region of integration from Fig.5 - 
& 3 +2, n=n ® 


em this transformation, it follows that 
AO'A' B AQ”A”B*® je 
Evidently 
0,0 > Ok, 
Ad, — Adı,o 
Bio, > Ba, 1) 


The JAcoBIAN of the transformation is 


> 
N 
—— 1 
Fr 4 
Therefore, the integral equation (46) becomes 
*ı 
h 2 ma? Ad [x &, ma (29, — &)] dE 
Gil, y) = — —— ı - „—ild ih 2 
En, = [ mean CM er 
Mı 
‘To write equation (48) in a more convenient form, let 
AR 49 
N = &a m , ( ) 


BET, Transformed region of integration from Fig. 6 


ER m 


=  eyuation (48) becomes 
# - * 2m Ale&,ma& 2m — 1] EdE, dr 
Gi, yı) = — ea Er ra ee 51). 
Era ef] ymab m VP Er 
Jefii 7; e &* by +7 
. gr u De a ae EEE (52) 
Then, equations 61 and (52) become 
Rn *”g FE ’ 
1 Ag [&, m &* —1 * 
2 er) = A DIE EI (53). 
\ ur £ | 
Be air ar | 
ee ö 6 
 Differentiating of equation (53) with respect to x establishes that 
E E> GEM __ 1 "Ar [z, mx (2, — 1)] dr, et 
BE i 02 - 2nmz Fri Zt BAER 
> - en er = Au 
E < we i 0 
ir, .  Alm22%—1)] = 2 BE DE ee en, (55). 
Br - Then, equations (54) En (55) lead to 
ET (56). 


ET - PP) di 
g= rn | ( 1 


‚Note {hat Ge, ) can can 


im 


Define 


A delta-wing surface of the form of a trtahle cz 


formation (58). It can be seen that the boundary lines, 2: =+n 
as shown in Figure (8). Further, if x = h, then x* = h. The p« 

is mapped by equation (58) into line «* = 0. Hence, 
surface of the delta-wing into a rectangle with the boundaries 


and 


Fig.8. Transformed wing Dar from the Fig. 3 


I 
Therefore, equation (54) can be written as = 
2 $ - RE 
9GH(x* OGHat, ı 1 u 7 d ; ö 
n=1 Se 2 ee 
2 2 2 
0 
Define 
. = 
2 2 d 
{n(y*) = ze u Mr 2 (60). 
k na D7 ") 
v 
Then, ' 
OGH(x*, y*) N An(x*) 
u 2 ee N. (61). 
Assume “ 
(#) 
A| — 
h nn ı* 
nr Fr BI Tr Re De (62), 
2m (w) 


where h is the distance from the origin to the base of the BR we have 
Gy (@*, y 


N 
Bra ı sin 


n=1 


Hence, from equation (63), it follows that 


h S i 
_ __2 f oGlar,y*) . nna* 
Ky*)n = — ı | Ge er (64). 


v0 
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Moreover, we have to solve the integral equation (60). The main difficulty in solving this 


_ singular integral equation lies in the fact that the integration has to be considered in the Han«a- 


MARD sense. Much effort has been used to solve this equation by reducing it to a lower order 
of singularity by means of integration by parts. However, the method of integration by parts 
implies vanishing of the function Y on the boundary. It was proved that by such a procedure 
the boundary conditions (1) = 0 are not satisfied. Hence, this shows that integration by parts 
is not legitimate. 


In order to solve this equation we use the principle of superposition. 
Define 


A a ER ee br (65) 
Then equation (60) can be written as 
“ R “ar 
v(Ö + N 
2 
+ = Er A Te 66 
Ky*) Wr 5 (66) 
it 
Define 
SR et C+1 M 
a ea D 
and 
BT ES Ad ee (68), 
where 
LU 2 43 u ae see (69), 
(with K as a constant whose value is to be calculated later), and /,(y) is chosen so that 
set yr BR En 
* — re ee ee er teT Ya . 
| [AU yı Fe y* 2 Y | 
Also define @.by 
yF = —cosp 1 (71). 
Then equations (68), (69), and (71) lead to 
f- eosp) =—Kceosp +f(l—C0SP) .. ».. 2.2.2... (72). 
However, from equation (72) it follows that 
re cos 9) cosp dp = —K [cos pdp. 
‘ 0) N) 
Hence, 
2 
Ks iite EOS CO OId Unna (73). 
n 
N) 


Moreover, from equations (66), (67), and (69), it follows that 


* +1 


le )ae 


9 


Kr N ee I ld 
nen 7 
and 
* 1] 
p, E de 
*) — EN SE NER TE 0 No 17 RER (75). 
f(y*) = 2 (y+* — 5): 
As a first step in the solution of (75), this equation is written as 
* + 1 
rl 
1 5 (p* ” & . ) + (76) 
5 hy") = (y* — 8 B 54 


e > Br 
Br E 
r Fun na Wal 
re .. ER 2 Fe Kari 


ie: Ru Ser 


A A Ar“ 


The solution of this equation is given by 


ne ar 


D 


E.. Er 
3 : Ey) N—® Ber 


R nm per 2m: 
4 —1 =, 
Define 
= y* = — cos . Eu 
y " a . (83). 
— Co 
Then from equations (70) and (82) follows that 7, is a constant whose value is to be BE 
lated later. 2 
Therefore, from equation (80) follows > 
| +1 EN Fur) yT- — pr? d£* BE 
n nn a I mg 4. 
Ze )- ae "2 ) Ken [ Ye 
—1ı 


2 


Then &= + 1 implies %, = 0; hence 


Ku, L=0: 


Therefore, equation (84) can be written 


u) = a | (| ar Vi gr as (85). 


(*—-y) 9 ee 


By the use of Cavony’s Theorem of principal parts, it can be easily proved that (85) is 
the solution of (76). 


a5 


| 


ET PETER EN 
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. After the integral equation (74) is solved, the superposition of this solution with that of (85) 
gives the complete solution of (66). Equation (74) can be solved directly, once the corresponding 
integration is carried out in the sense of HaDAamArD. Hence, equation (74) can be written 


Ru 


(E41) 0m) 


Ye 


IN { dd 
— y* = lim d ay* on 
27 Tun J ed) a ee We 
el 1 
1 
Y ( = Zr) de 
+ d NE ——— 2. ; 
(y* — a) (y* — 9° =) 
t y*-+e y*+e 
Hence, 
Hu 2 nz 1 
K + N 2er) 
39 re Wr pr ec 6 (87). 
Note that £ means that the CAvck#y prineipal part of integral is meant. 
Define 
= Wil 
ae EA a a nn (88). 
b= Y,(0) 
Then 
: 1 
ee sen 
ae m | 89), 
where 
5b a 


Hence, 
et a gr 
ie 1 er F(y*) Y1— y* 
Y en Ti er a 91 
( : ) ia" | a = 
—1 
where 
-1 
, DEN 
ja Be) EN N (92) 


After substitution of equation (91) into equation (92), it follows that /\ is an undetermined 
constant. Therefore, from equation (91), it follows that 


E tl —— 
FENG .D, de* 1 F(y*) yı a br nd, 
2 Be) en | Yen nr = ( Er — y* ") Y1—c*: ee (92): 
| —1 —ı 


2 —1 2n—1 +1 NE 
na a | Trwrie ,.) a h 
Fl) = — | Vi gr 7 | ( Er yi dy ir +G .(B). 
en] —1 —1 


Hence, 


702 


I (95) and 1.06) lead to 


u 


Yu) = A 5 Yes F 


a—cos!(2n—1)/ 


v 
a—cos!(27,—1) x 


a dr 
a Tu ya + c0s p) (cosp — cosy) 
0 


Therefore, 


za 


dy—+C, URÖIE AL u r (97 


a—cosH(2,—1)/ a 5 : ee 
coS og Sın { . 
Yun >; [ IT ap 


0 
a — cos-!(27,—1) 


n (1 + cos p) dp 


cos y — cos p 
En.“ 


ae Te nee 


m cos y — C0S@ Mi 
0) 0 
7 — cos! (2n.—1) Ei; (i ) d 
a a A 
en | Sen IT a ren 
v N) 
However, 
RR 1 1 
cos p sin = 76059 —7005pc0529 else ee 
and 
cospcos2p = 2 (cosp+cos3p) ..... 
Then, equations (99) and (100) lead to 
| cospsin’p = ei cos p — & cos 3 
[17 haar a ziie BE 
The use of 
_cosnp 5 _. _asinny (n=0,1,2,...) 
cosp—cosp 7 sin y O<y<a 
in (98) leads to 
Fl) = +2 — [ar 0522, — 1] — Ey — Im — cos! n— 1] 


2 — 005-142, —1) 


sın Y It 


K sin 3 b 
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Moreover, 
sin3y A 
ip = 2cos?y + cos2y. 

Then, 

nr — 605 1(2,—1) in3 
sin 3 y K a Br: 5. 
37 en dy = Ep [r — cos! 2 — 1] + 37 sin2 [a — cos!(2,»—1)] (103). 
0 
Hence, 


IF = Re 
Fo) = [ra — cos! 2m —-— DIT, +b—- + z., un 2 [r — cos!(2,—1)] + C, (104). 


From equations (88) and (104), it follows that 


NIC 
ER = a ige: (105). 
Therefore, j 
I, e er he 
Ale De [r — cos (27, — 1)] — 3 sin 2 [cos!2 (m; —1)]+Cı - . . (106) 


is the solution of equation (74). Since /, and C, are not yet determined constants, equation (106) 
must be proved to be the solution of equation (74). For this purpose, equation (74) can be written 
in the form 


”1 
F) dn 
*— el ET ee ee N 107). 
| + nm — DR Er 
(0) 
Equation (106) and (107) lead to 
K ae dng I Een (2n—1)dn 
4 Droemer = y* —- 2% — )P 
Be sim 2. IE @n— 1] dns 
8 ie RT RI VE WERE Od (108). 
a) 


Note that each of these integrals has to be evaluated in the HADAMARD sense. 
Define 


“1 “+1 
_ ff Benni: til FIR 
= 2 5 | Maar es) air eao)N Kr = Bu rs 
0 —1 
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M, = & ee (110), 
—1 
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Ten) Wr em—Dr ee en 
, —1 
Therefore, 
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M, = a lim 
wre .)0 nt „e 


or 
T, 1 
—-I— Te a a 112). 


to see Er equation 1a) leads to 
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Then, by integrating (117) by parts it follows that BERN en 


en! K ee. 
5 en yr—Hyi—a 


wi Te 


Or 
RE: rn est hr au x FASER 
In. (—H)y—a 5 j 
Evidently, 
x En ar RE 
Se 4m w—d Ke an 
or 
K f costy dy | er 
M, = 1 = cos y — COS BB I Eee Zn er a: 3 
Hence, | S 
M=—7 009 ER (122). 
Therefore, I . 
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Equations (107), (112), (116), and (123) show that 


G, = 
er (124) 
Hence, equation (106) becomes 
Kerr 
vn) =— gu RS et) ae (125) 
or 
K Be Ya ann 
Pl) = 7,12) er DD Se (126). 


Therefore, equations (67), (85) and (126) lead to: 


Disker, 
mg 


N +1 EDER 
K Er A Yylyr de* 
”(n) = 12 1-29) —-@,— 1% + Im { in on dy* on (127). 


ua 


This solution must be expressed in the forın of (x,, yı) coordinates. Equation (55) can be written 
as 


n=1 


Zn! eu 
EN 1 (y*)n Y1 — y*? de 
Ai man N) 3, Ya 15 | AT BF: erg 


| al 


K 
yon -HDN-OR-—T Zerl28): 


Equations (52), (49), (48), (44), and (128) lead to: 


Yı 
Mı &ı 


+1 Eu Rn Turn 

1. 1 *) /1 — y*2 de* ne 

A = j a 2 a no Yı Fr z| z ae) 
n=1 2 — a 


—1 —1 


(129). 


Equation (129) is the basis for the general solution for transonie case, since in order to 
evaluate Af(z, yı), ..., we must solve an integral equation of the same type as that given 
by equation (16). 


Therefore, the next step is to evaluate equation (17). Define 


1 Ayla, Yı) dx, dy, 

xy) = — - ERTL TEE 0). 
EA (130) 

D, 

Then, the equation (17) becomes 
1 A% (2, y,) da, dyı : 
uU) = — — = : RENTEN: 151). 
SE 1 Il eh) 
D 


Equation (131) is of the same type as equation (16); therefore, the solution of (131) is given by 


N | 1 ma Er (y*) UL 2 de* R ie - | 
1 Gym —yY € a Yı Yı 
*/n he FREE A wE * ee = & 
Aslar yı) = 27 2, Anlcı) jE N j Er y* dy ı-cr 2 e 2 Vi E a | 
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where 
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and 


R=—, [atesnenspdn u a Re 
u . 
v0 


Moreover, knowing Af(x,, y,) and A3(x,y) we can easily determine Ara 7,). Ine 
process of evaluation of A}(x,, y,) can be extended as far as we please. ‚However, it is necessary 
to point out that the terms in the series expansion of the kernel, equation (11), decrease rapidly 
since 2 — 2, > ß(y—yı). From the point of view of accuraey, two terms in the series expansion 
of A,(x,, yı) will be sufficient. ar 


Hence, equation (18) can be written in general form, namely 


—1ı 


ae Imre 
1 «| Ne nt (yrHyn +) —y** ar 
A, 40) = ef 2 l >, PT? Anl) N j u L : &* u dy* yı — + 
j=Ul In= En 


Equation (135) represents the coefficient of the linear term in the expansion of A of a source 
function for the transonic case for the non-steady three-dimensional flow moving around a delta- 
wing. The solution, for the value of ß greater than zero, can be obtained using the same mapping 
technique and procedures as were used in the solution of (19) and (20). 


CGonelusion 


A systematic presentation of the solution of the integral equation of the non-steady 
motion of the delta wing in supersonic flight has here been done. The new method of solution 
of this equation consists of the perturbation of the domain of integration. The solution is given 
for transient condition, but no difficulty arises for obtaining any other supersonic solution. 
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The Stress-Strain-Relation of Elasticity and the Measure of Strain 
By M. Reiner *) 


Die Spannungs-Dejormations- Beziehung der Elastizität, welche vom Deformationsmaß unabhängig ist, 
enthält den Deformations-Tensor in nullter, erster und zweiter Ordnung und eine unendliche Zahl von skalaren 
Parametern. Wenn der Verschiebungs-Gradient infinitesimal ist, lautet diese Beziehung folgendermaßen 
$j en (AI+5I) dr +2 mein +4 Me Ein en mit vier konstanten skalaren Material-Parametern A, ö, u 
und We- 


Wenn die Parameter ö und u. verschwinden, wie dies in der klassischen Elastizitätstheorie angenommen 
wird, dann kann das elastische Verhalten eines Materiales nur bei Einführung eines besonderen Defor- 
mations-Maßes e;r adäquat beschrieben werden. Dieses Maß kann als das ‚‚natürliche‘“ Maß des Materials 
betrachtet werden. 


The stress-strain relation of elasticity which is independent of the measure of strain contains the strain- 
tensor in zero, fürst and second order and an infinite number of scalar parameters. When the displacement 
gradient is infinitesimal, the relation is reduced to 5, = (AI+ 6TI) 5a +2 u en + Aue Ein Eur with the 
four constant scalar material parameters }, 6, u and ie. 


When, as in classical elasticity, the parameters ö and u. are assumed to vanish, the elastic behaviour 
of some material can only be adequately described by introducing a specific measure of strain eir, which can 
be considered as the “'natural’’ measure of the material. 


CoOOTHOIIEHHE NeBoPMauHuH — HAIPSGKEHHA B YUPYTOCTH (9JIACTHYHOCTU), KOTOPOE He3aBH- 
CHMO OT MepbI NePopMAalmH, CONeP5KUT TEH30P — Nehop;kammnu B HylleBoü, IIepBoü, BTOPoH 
CTeMEeHH MH ÖECKOHEYHOE KOJIMYECTBO CKRANAPHbIX HapaMmeTpoB. Ecam TpammeHT TTrepeMemeHun 
ÖecKOHEYHO MaJl, COOTHOIIEHHE YMEHbILEHO R 5 — (Al + ÖL) dr + 2u&ir + A ucEia Ex C YETBIPBMA 
HOCTOAHHBIMH CKAJIAPHbIMH IAPAaMeTpaMmu MaTepnada A, Ö, u, He. 


EcJH Kar B KJIACCHYecKoON TeOPUH YUPYTOCTH NOJIATAW TUTO Ö U u. HCUe3HYT, TOTAA yupyroe 
IIOBeIeHNe KAKOTO-AIM6O MAaTepmasla TIpaBuJIbHO ONMCAHO, TOAIbKO ECJIU BBOAAT OCOÖEHHYIO 
Mepy NehopMaumm. JTa Mepa MO}KHO CYUHTATBb RAR ‚‚ECTECTBEHHYIO‘‘ Mepy MaTepmalla. 


1. The stress-strain relation of classical elasticity is expressed by the linear equation 
Sr =Aldir + 2u 6% ae Er RN 


where s;; is the tensor of stress, e;; the dimensionless tensor of strain, I its first invariant and } 
and u are two constant scalar material parameters. The stress can be unequivocally derived from 
the concepts “force” and “area” as defined in Mechanics and Geometry. The concept strain is 
derived in classical “infinitesimal” elastieity from the concept “displacement” (1;) as follows: 


It is obvious that a homogeneous displacement does not involve strain. There must be 
a displacement gradient /}; = u; ;. However this statement cannot be reversed. There can 
be a displacement gradient in the absence of strain. This is so in the case of a rigid rotation. 
If the rotation is infinitesimal, it is defined by the tensor, 


1 
Vik = 5 (U,r — U; ;) . . . . . . . . . . . . . . (B2r. 


Therefore u; — (li, — U;,;) can be considered as defining an infinitesimal strain 
| (1.3) 
&r > 5y (u;,x — m Uri) BE Fe" en : 


the gradient y;r = U; beeing infinitesimal. 

These considerations cannot be applied when the strain is “finite”” because (1.2) does not 
define a finite rotation. Here the procedure is as follows: 

Consider a spherical part of unit-radius of the body under consideration. This sphere may 
be finite when the strain is homogeneous; when the strain is not homogeneous the radius of the 


sphere must be infinitesimally small. Through straining, the sphere is converted into an ellipsoid 
with half axes 1, I, I,, or in general !. Any function of I which vanishes for 1 — 1 can be considered 


*) Israel Institute of Technology. 
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as a measure of strain. The following measures have actually been proposed by CAucHY, GREEN, 
Hi&ncKY, ALMANSI and SWAINGER!): 
ee =1—1, 


ee 3 nz 
EN NAzl1, 
ı 
dl = er 
= (Zeini=dD a Br A a 


1 


1 


=. 


ed 


ig) -gum4 


= 1 ddl 4ır 


In order to appreciate the differences between these measures, consider a rod extended to 

twice its length. Then the different measures give 
e — 100% , es = 150% , eH = 66% , ei = 375%, ee = 50%. 

Therefore, if we keep Equ. (1.1) as our stress-strain relation, we should find different 
stresses for the same dimensionless strain, depending upon which measure of strain we use. This 
is a paradoxical situation. A worker in this field has said on occasion: “(this) cannot ‚possible 
hold... (To say) that the stress... depends on what tensor you use to measure strain.... 18 
wrong. The formulae may indeed look different in different formal structures, but the result, 
the numerical values of the stresses, must be independent of the kind of tensor, the theorist 
chooses.” This sets a problem with which we shall presently deal. 

2. If we take a rod of some specific material, say rubber, extend it by 1% and determine 
the stress necessary to effect this extension, then extend it by 100% and determine again the 
stress; the numerical values of these stresses will of course have objective existence, independent 
of the measure of strain chosen. We can easily derive from Equ. (1.1) the reverse strain-stress 
relation. 

j m Sik >) 
er >= Fu@A +3 u) "ed "Zu u. ne ee (2.1). 

Introducting the values for the stresses as found experimentally, we shall find values for &; 
which will in general not be 1% and 100% in which case CaucHys measure is applicable, but 
something different and the question is whether there is among the measures listed in (1.4) any 
one which gives these results and therefore applies. 

Let us assume that the stress corresponding to an extension of 100%, was 66 times the one 
corresponding to an extension of 1%, this would show that the measure to be used may possibly 
be the logarithmic or the HexckyY measure. But if the material is rubber this will not be the 
case. War (1942) has proposed a different measure claiming that it is applicable for rubber. 
It is 

er =- (RR —Iid)=2e +0, 


Therefore, if we start from a linear stress-strain relation with constant parameters, such 
as (1.1), we are not at liberty to chose our measure of strain. There may be for every material, 
in accordance with its actual behaviour, a different measure appropriate for it, or what we may 
call its “natural” measure. Such measure can possibly be found, as has been indicated above, 
by experiment. It can possibly also be deduced from considerations of the internal structure of 
the material. WALL arrived at his measure in this manner. Nevertheless, it must be admitted 
that this situation is not satisfactory. The ideal would be to possess a stress-strain relation which 
would be “independent of the kind of tensor (of strain) the theorist chooses”. 

3. Reiner (1948) has solved the paradoxical situation as follows: 

If we relate the tensor of stress in a general manner to the tensor or strain by 


Ss: = Fan) a Be We Ne 


and assume that the material is isotropic, then, because of a general principle of dimensions, the 


funetion F must be expressible as a sum of symmetric tensors of second rank, multiplied by 
scalars, thus 


\ == cr R 2 Rp r p 
Br Si = do dr + Bı 5 + Ban Eur + Or np Feen (3.2), 


!) For reference see TRUESDELL (1952). 


a» Ei r 
0 „ \ 
MET 


a a a in ln aan ann u, 


Lu 5 


EN 


SU 


where the 5 are in general functions of the prineipal invariants I, II, m of er. 


j 
« 


in Relation of Elastieity and the Measure of Strain 


” 


Drawing upon the CayıevY-HAmıLroN equation of matrices, he has shown that Equ. (3.2) 
can be reduced to 


SE = (Foı IH Fol On er + Friae Aigaoe Fo + 83), 


Pr 


where FuF, and F, are functions of the three principal invariants. Following Tnunspauz these 


can be expressed i in power developments 


Fa = Faooo + Frasoo I + Fnzool? + Faoio II + Frag I + Frau DIT + Frog IT ++. (3.4) 


so that (3.3) becomes 


Sik Rs Id; + F'o00 ex] + [(Fozoo I? + Foo ID) 5x + Firoo I er + Foo Ei Car] 
+ IFosoo I + Foo III + Fooo III) 5 + (Fi200 ? + Fioi ID) &x + Faro l&iu rl +: (3.3). 


Where the expressions within [ ] brackets are of the first, second, third and higher orders 
respectively. 

Note that the term Foo must vanish as the stress vanishes with the strain. 

We see that the stress-strain relation (1.1) is a degeneracy of the more general equation (3.5) 
in that it assumes the strain tensor in the zero and first order, while the complete equation 
requires that it should appear in the second order also. Moreover it assumes that there are only 
two scalar material parameters, both constants, while in the complete equation (3.5) there appear 
theoretically an infinite number of scalar material parameters. Under such conditions Equ. (3.5) 
will be adequate to express any stress-results in any experiment using any measure of strain, 
with, the only difference that the material parameters Fyop will have different values for different 
measures used — arestriction which is acceptable. It may be mentioned that in the early attempts 
to describe rubbery elasticity, while keeping Equ. (1.1) and CAaucays strain definition, a variability 
of the Young modulus was assumed which meant introducting parameters Fynop- This procedure 
was ultimately not successful and had to be given up. 


Thus a complete mathematical solution of the problem is attained — but it will not be 
satisfactory to the physicist. Many simple relations involving a small number of parameters can 
be expressed by infinite power series with an infinite number of constants. This is not what the 
physicist is after. As said before, the parameters will have different values depending upon the 
measure of strain used. Some or many will vanish. The physicist is thereby thrown back to the 
attempt of finding such a measure, the natural measure of the material, that most para- 
meters vanish, except a minimal number of them. 

There are several such attempts in the field of finite rubbery elasticity, in addition to the 
one by WaArr, mentioned before. No such attempt has been made in the field of metals, where 
the strain is infinitesimal. 


4. We must therefore put the question: what are conditions when the strain is infinitesimal ? 
In this case higher orders of 1— 1 in Equ. (1.4) can be neglected against 7— 1 and it seems that 
all strainmeasures are reduced to the one CAucHY-measure which is also embodied in Equ. (1.3). 
This, however, is correct only, if we consider the principal values of the tensors as written 
down in Equ. (1.4). The contention that the difference between the different measures vanish 
in the infinitesimal case, is not correct when we consider for instance shear involving rotation. 
To show this, we must draw upon the complete tensor expressions. 

The tensor expressions for the two measures of GREEN (G) and Armassı (A) have been 
derived by MURNAGHAM (1937). 

They are 

2; = ul’ + iR + @1) (a): ] 


„t4l). 
a 7 (P, ai) (Fax) [2% 


In these equations an index placed to the left refers to the “initial”, one to the right to the 


“final” state, so that 
vu —ı 
An. N N 
T 1 
wh = ou;/drX ; T’;r u o1;/0X% f 
Consider simple shear given for instance by 


eg: Be me... nn ka). 
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Therefore 


Ge 
Ta=|0 0001 “AP... Sa 
000 


The straintensor is in the two measures 


r Rab r|® 1 | Fi: 
me=— 1027-000 ; Ah=z |! —Il aa a SET PN? 
000% io oo) 


The stress must be refered to the final state which is the state of equilibrium. The ÄLMANSI- 


measure refers to this state, but when the GREEN measure, which refers to the initial state, is 
used, we must carry out a rotation which brings it into the final state. 


The rotation matrix is . 
'". Cosa’ Bin 


Rır | . (4.6), 
0 0 1 
where (compare Lovx, Art. 3.) 
tan Tan II EDEN ER 
This makes 
r Tau 9 
ee nn 
- 119..0..0 
The three invariants I, II, III of these two measures are 
ee Tre r:2, 
ae en 
Ir = IM — IM, 1119 = IF =0 


Now assume J’to be infinitesimal when we write it y, then y? is infinitesimal of second order. 
Nevertheless it is not negligible in either 


ee a ee 


nor 
A ee Re 


For a quantity to be negligible it must be an inerement to a quantity of lower order and 
in addition its magnitude must be less than the error in the measurements of the quantity of 
lower order. For instance, the maximum elastic extension of a metal wire is of the order 10% 
and therefore “infinitesimal” in comparison with unity. The maximum positive and negative 

“eross-strains’’ of Equ. (4.10) and (4.11) will therefore be of the order of 10%. Nevertheless these 
are not negligible, appearing in their own right and not as increments to first order terms as is 
the case with those of Equ. (1.4). If simple shear is applied to the case of torsion of thin wires, 
they become manifest in extension, transversal contractions and cubical dilatations, as was 
observed by Poyxtına (1909, 1912). 


9. REINER and Rıyteu (1958) have applied these considerations in order to find the form 
which equations (3.5) will take in the case of infinitesimal strain which we denote bye. 


When the first-order term of (3.5) does not vanish, all other expressions can be neglected 
and we have 


Sy = Fam 10 + F 000 &ir ee (DER 
We consider two such cases: 


(i) We first subject the specimen to an isotropie pressure (—p) so that 
Ojk pP Ösk nf de fi 00 uk 2 er ER (5.2). 
This is connected with a cubical dilatation of isotropie displacement gradient y so that 


I=-3y= 5, 11-372, Te ia ae 
and (5.1) becomes v } Y (5.3) 


P= 3% Foo + Y Fiooo = &r (Faro + El) en (5.4). 


ee) 


a et ne 


tee 
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We find accordingly that 


KR a 2 Br a De ER RER Tr Be 1 ee 7 (5.5) 
with x the bulk modulus. 


(ü) The second experiment is simple pull when for instance 


000 
Si BO a ee (5.6) 
006 


Ba Vs a Re (5.7) 
and two transversal contractions 
& nn €o = Be ET ee ee ee en ar (3.8) 
with 
I=py—2y =8&; 1l=y(—2y); INT arten (3:9): 


Equ. (5.1) now yields 


ee ES N IT (5.10) 
0 = Foo (yı — 2 Y1) + Fiooo Yı 
Introdueing Youngs modulus E and Poıssoxs ratio » so that 
or =E; N ee (5.11) 
(5.10) bescomes 
Fan(! —2r) —- Em’ =0, ! 5.12) 
Foo(l — 29) + Fi = E 


from which 
Fool+»)=E, 
Ds an es (5.13) 
Foo = #— (2uj3) = A 

as in classical elasticity. 


(iii) Conditions are very different in simple shear as treated above in $4. We have in general 


Erz Ey 


Eyz Eyy | 
Therefore | 
| 22 + &y Ery(Exz Pr &yy) 515 
WE A Re ee) 
Ery(Eas + Eyy) Eyy I En 
If &,, = yl2 is of the first order, e,, and e,, are of second order and therefore 
| y 2 
| 9 Eor F Eyy | 
Eia Eak — = | | . (9.16). 
= R 
(ee SE Eyy D7 | 
We have now 
ar ER yet BE 47 
Out Fuel + Ems lee ta) Fr ne... ON. 
On appropriate approximation this yields, considering that I is of second order, 
ei . (5.18) 


asin classical elasticity. 
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It is different with the normal components. We find from (3.5) taking into consideration 
that &,7, &y, and I are of second order and neglecting third and higher order terms 


« F. y? F. ” 
S=ill+r2u, oo + 2000 4 


y? z X is 
Sy=Al+ 2uey— Foo g + Fa >» er Pe Be RE une Gr 


5, AT — Fo 7 


In general, in order to produce simple shear it will therefore be necessary to exert, in addition 
to the classical shearing stress s,, of Equ. (5.18), an isotropie “cross”-pressure of the amount of 


2 


Po=Fg—Aal 2.22. 2 6) 


% 


When such cross-pressure is absent, the volume will expand. This phenomenon was des- 
cribed by Revyvoros as “dilatancy”. We may therefore term the factor Fy, modulus of dila- 
tancy d. This should not be confounded with KRONEckEss delta. 


In addition, there will be normal “cross-stresses”” $7,., and Syye 


RT a u EI 
Syvo = 2 U Eyy tr Hey” 


Where 
Fzo00 = A He E 


and „.isa modulus of cross-elasticity. 


6. From the foregoing it appears that the general law of infinitesimal elasticity can be 
written as 


Se =(AI+5 MD) dr + 2uEa + LUeEis Ear u (6.1). 


Where A, ö, a and u. are four constant scalar material parameters. First and second order terms 
only are present, the second order terms only in the absence of first order terms. Note that in 
general I will be of second order in simple shear. Third and higher order terms are negligible. 


Any measure of strain can be used for eg, but the numerical values for the four material 
parameters will depend upon the measure chosen. The attempt to find a measure for which ö 
or,and u. vanish is still legitimate. It should, however, be pointed out that in the absence of ö 
and w, any measure for e will require the PoyxTisg effect to be an extension and not a shor- 
tening. 

Therefore, should a metal wire in torsion show a decrease in length, the complete equation 
(6.1) must be applied.) 
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°) PoyNntingG tested two steel and ‚one copper wire. Both showed lengthening. Foux and REınEr (await- 
ing publication) have tested 11 steel wires. Of these 10 showed lengthening and only one shortening which 
could be attributed to an inhomogenity of the material. i 
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ıktion von Dirac in der Wärmeleitungs- 


Der Hauptvorzug, welcher zu einer immer breiteren 


Anwendung der ö-Funktion von Dirac in der Analyse 
von klassischen Problemen der Physik führt, ist die 
Einfachheit und die Schnelligkeit, mit der man die 


Endresultate erhält. Die ö-Funktion ist eine Analogie 


_ des Symbols d,,,n von KRONECKER und stellt seine Ver- 
 allgemeinerung dar in der Hinsicht, daß sie die Ortho- 


normalität des kontinuierlichen, während d,,,n nur die 


 Orthonormalität des diskreten Spektrums bedingt. 
In dem Artikel wird die Anwendung der ö-Funktion 


auf die Theorie der Leitfähigkeit in einem drei- 
dimensionalen nichtbegrenzten und begrenzten Raum 


angegeben. Indem man die Analogie der Symbole 


ausnützt, wird die Lösung des Problems der Leitfähig- 
keit im begrenzten Raum angegeben, welche der Lö- 
sung von PATERSoN [1] für einen dreidimensionalen 
unbegrenzten Raum ähnlich ist. 


a) Der unbegrenzte Raum 


Die Untersuchung zur Lösung des Problems der 
Leitfähigkeit im unbegrenzten Raum: 


20 


ar rotrreEyne) Ha MD, 


O0722,0) 9, Ze (2), 


kann man in eine Untersuchung der Aufgabe eines 
einfacheren Problems überführen: 


ru — 5,8) 2 0.8), 
u, Y»%, 0) —- UylR; Y; x) gen. on (4), 
durch Substitution von 
t 
u= 0 exp I S ve) a) H 
> Fila: 


t 
g=ypexp = S vw) a) 
Ö 
Schreibt man die Lösung der Gleichung (3) mit der 
Anfangsbedingung (4) wie folgt: 


u(z, y, 2,8) = //[SS (&,y21;%0,Y', 2,t) g(@', y’,2’,t‘) 
Badahduykdzädte re ee (6), 


so erhalten wir — mit Hilfe der Formalität der ö- 


Funktion — den singulären Teil der GREENschen 
Funktion G, aus der Gleichung 
G4=— Lö —a)öy—y)d@—r) (M 
oder in entwickelter Form 
ib 

A= Omi 

X SISS In ei h@—r)Hky—v)Hks@— 2) HiNHE) 

Es) a u ML Er (8), 


wo. nach (3) 


(0) 
ner 


22 2 8 
= LT: 


0x? 5% öy? 


Somit ergibt eine umgekehrte Operation in (8): 
In aan) + kW y)+ ka @— IHRE) 
elkı@—r)+ kay—y')+ ka @—z')] + Ehlt—t) 


= 10), 
x ehe tih m 


wobei 
P=M+M+R. 


a 
itteilungen 


KLEINE MITTEILUNGEN 


Be 
Zr Rn 

er re 

IE 


» 


- Wenn wir (10) in (8) einsetzen und über % inte- 
grieren, erhalten wir 


u 
om) | 
x [Sf er rUudH Man) + Lu N)+ he) 
Un dk dh,- ou a, 1), 


und nach nochmaliger Integration 


LRZEEUZEUNERZEZI 


Rear rn 4#(- 1) 
= (2x Ya — r)® fürt>t’ 
Den für t<t 
(12). 


Die Funktion G, entspricht der Temperaturvertei- 
lung einer Einheitsquelle, die zu der früheren Zeit 
i — ti’ wirkt im Moment t, was einer Wirkung mit Ver- 
zögerung entspricht. Dieser Charakter der Wirkung 
geht aus der Form des Operators Z hervor, für den es 
keine Lösungen mit beschleunigter Wirkung gibt. 

Wenn wir (12) in (6) einsetzen, finden wir die Lö- 
sung (3), welche dem singulären Teil der GreENschen 
Funktion entspricht: 


ih 


us(&,Y, 2,1) = 1 | wa 
3 (2% Yr) ) Ya— 
{0 


az + We) 
I 


er en l)dehdyadzeidtie sn. (13), 


welche die ‚Wärmestrahlung‘‘ von Quellen der Inten- 
sität g(z, y, 2, t) bestimmt. 

Die Integrationsgrenzen gehen sinngemäß aus den 
Bedingungen (8) hervor, welche das Vorhandensein 
der Funktion @G, bestimmen. 


Für eine Quelle, die im Moment t’ = 0 wirkt, stellt 
die Funktion (12) die Elementarlösung dar, deren 
Superposition erlaubt, Lösungen aufzubauen, welche 
die Rand- und Anfangsbedingungen erfüllen. Indem 
wir diese Eigenschaft benützen, kann man eine Lö- 
sung ermitteln, welche die Rand- und Anfangsbe- 
dingungen erfüllt. In unserem Fall kann eine Lösung 
gefunden werden, welche dem nichtsingulären Teil der 
Greenschen Funktion G, entspricht und die homo- 
gene Gleichung LG, = erfüllt. Im Falle der Be- 
dingung (4) haben wir 


1 
(2% Yn D)° 


» RT RVM > 2) 
= fee, y',2’)e ie 


x da’ ER RE (14). 
Auf diese Weise stellt die Funktion 


una, = 


u2,y, 2,1) =, 20) +u,(,y2%t (15) 


eine Lösung der Aufgabenstellung (3), (4) dar. 


Die Einfachheit der Methode zur Lösung der ge- 
gebenen Aufgabe mit Hilfe der ö-Funktion ist leicht 
erkennbar, und zwar durch den Vergleich mit den 
Lösungsmethoden der klassischen Analysis. 


b) Der begrenzte Raum 


Um in diesem Falle einen singulären Teil der 
Grernschen Funktion zu erhalten, drücken wir die 


ZI 
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ö-Funktion in Form einer genormten ‚Reihe der 
Funktionen wn(x, y, 2) aus, welche die Gleichungen 


L wm, 9,2) = —A um, y2) - - (16) 
erfüllen, die Bedingungen der Orthonormalität 
b 
{5 From, 4,2) unla, y, 2) da dy de — Ömn (IM), 
000 - 
sowie die Randbedingungen, gegeben durch die Be- 
dingungen der Aufgabe für die Funktion w(z, y, 2, t). 


Das dreidimensionale KRONEcKER-Symbol nimmt 
dabei die Werte an: 


Ömn = Öm,n, Om;n, Öms n; 
1 form en, men,m=n,t (18). 
“ 0 sonst für andere Fälle 


Das bedeutet, daß wir die GREENsche Funktion mit 
Hilfe der Gleichung (7) suchen, in welcher gilt !) 


Öle — a’) öy — y') de — 2’) 


= II I uml& y, 2) Wmle’,y',2’) (19), 
MM; Ms 


ut) = RR dk... (20). 


Wir haben also 


1 Wmlz,Y; 2) Wm(x’, were 
g 2x ö 2.2 dk 
e nn Mxt+ik 
= BIRNEN Wmlz, Y, z) wmlx, y', 2’) As »2(t—t') 
Mm, M; My; 
era le. 
und daraus 
abet 
u,(r, Y,2, 2) == >> f SS wntz, Y; z) Wmlx', y, 2’) 
Mm, MM O00O0OO 


nnd AL) ur, y',z',t') dx’ dy’dz' dt’ (22). 


Die letzte Gleichung löst das Problem der ‚‚Wärme- 
strahlung‘“. 

Um die „Wärmeausbreitung‘‘, welche dem nicht- 
singulären Teil der GREENschen Funktion entspricht, 
zu finden, muß man die anfängliche Bedingung (4) 
durch die momentane Quelle der Intensität 


90(%, Y, 2,1) = uslx, Y, 2) Ölt) BEER 


!) Wenn manim eindimensionalen Falle die Funktion /(x) in ortho- 
gonale Funktionen wm) zerlegen kann, dann ist im bezeichneten 


Raum f(x) = FA LER GC) PR wobei 
Mg 


= F 

1 fürmen 
w, (x) w,(r) de = ö = 

{ TE mn R fürm#n 


sowie 


a 
An = { i(&) wm) dı 
und schließlich 


a 
I) nr Lem) Wr) de’. 


Diese Entwicklung kann man auch wie folgt ausdrücken: 


M a 
Ka) = lim FE Sfa)w,,(e) w,,(z) de 
Mo m, 0 se sn 
oder 
a 
a) = lim SI@) dla — a’) de’ 
M—oo (0 > 
wobel 


M 
dsya—ır)= LE we) we). 
m, 


Für M— oo wird 


ödr—r)= Lu (e)w_(). 
m, m m 


Auf diese Weise wird die Darstellung der ö-Funktion nach DikAc in 
Form einerorthonormalen Reihe (19) an Stelle der üblich angenomme- 
nen Bezeichnung nach Formel (20) begründet, 


Kleine Mitteilungen 


die Temperatur Null besitzen. In diesem Fall kann 


ei Ta a Werte ul2, 9,2) ge 
wähnleistet?). ar o % 
Diese Ersetzung erlaubt es uns, di (22) um 
ine Lösung, der anfänglichen Bedingung (4) ent- 
ee = Fan änelichen Bedingung (4) ut ’ 
erhalten wir nämlich:- 


abe E 
le, y, 2,1) -ZEr SS Sunte Y. 2) um(, y, 2) 


x" uular, y', 2’) de’ dy’de’ (24). 
Auf diese Weise ergibt die Summe der Funktion (22) 
und (24) die Lösung des Problems der itfä 
(3), (4) für den begrenzten Raum und für Quellen 
9x, y,2,t), die vom Zeitpunkt =(0 an wirken. 
Für Quellen, die früher wirken, beschreibt die Summe 
(23) und (24) die Erscheinung der Strahlung von der 
kontinuierlichen Verteilung g(z, y,z,f) und der Mo- 
mentanquelle (23). Daher soll man die untere Zeit- 
grenze des Integrals in der Formel (22) bis — © er- 
weitern). 

Zur Erläuterung der Symbole in den letzten For- 
meln kann man die Wärmeausbreitung im Parallel- 
epiped betrachten, dessen Kanten mit den Koordi- 
natenachsen zusammenfallen und dessen Rand- 


flächen 
=; y di WE 


man dem Ausdruck w„(z, y, z), der die Bedingungen 
von (17) erfüllt, folgende Form geben: 


/ nz: 
wm(x; Y, z) = J = zur = sın — 


ae a b o& 
wobei }? den Wert 
h mi, mi ‚mi 
2-2 (m +p+=) 
annimmt. 

Wenn wir die Transformation von (5) auf die 
Summe der Funktionen u, + v, ausnützen, kann 
man schließlich die Lösung des Problems der Leit- 
fähigkeit (1) und (2) für den begrenzten Raum in 
folgende Form bringen. 


t 
[ wer) " 


ö 


O(z, y,2,i) =exp 


x > >> > 9x, y’, 2’) wm(X, Y 2) wm(x’, y’, 2’) 


m, m, I, 


x et ger dy' dz’ 
t 


t 
z — S plt’)at' 
n e 0 
exp \ yet) “| ' a Ah 
ö Yet —t) 


) 0 


x | J >, 223 Wm(T, Y, 2) wm(x’, y’, 2’) 
ER a an 


Ele Opa, y',2’,t)de' dy’ dz’ dt’ (25). 
Diese Lösung ist der von PATERSoN für den Fall des 
unbegrenzten Raumes analog. 


Literatur 
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») Es wird angenommen, daß die Intensität g(x, y, z,t) also auch 
Ile Y,2, . For Tree mit x*/k die Dimension einer Intensität be- 
sitzt, wobei wie üblich x? und & Koeffizienten von Temperatur- u 

Wärmeleitfähigkeit darstellen. : = 


_ verfahren 


Zur Konvergenz der Differenzenschema- 


- 


1. Einleitung 


Vor einiger Zeit hat D. GAter [1] unter Benutzung 
einer independenten Fehlerabschätzung gezeigt, daß 
bei der numerischen Integration der Differential- 
gleichung y’ = f(x, y) die mit dem Apanms-Verfahren 
gewonnenen Nährungen bei kleiner werdender Schritt- 
weite gleichmäßig gegen die exakten Lösungswerte 
konvergieren. Im folgenden wird gezeigt, wie man 


für ‚‚Differenzenschemaverfahren‘“ zur numerischen. 


Integration von gewöhnlichen Differentialgleichungen 
beliebiger Ordnung diese Konvergenzeigenschaft nach- 
weisen kann. Dabei werden nur die aus den Fehler- 
gleichungen unmittelbar folgenden rekursiven Fehler- 
abschätzungen benutzt. Fragen der numerischen 
Stabilität werden nicht berücksichtigt. Um uns kurz 
fassen zu können, schließen wir — ebenso wie D. GATER 
— an die Darstellung von L. Coruarz [2] S. 111—130, 
an. 


2. Voraussetzungen 
Es liege das Anfangswertproblem 


ym = fla,y,y,...,y"»), 
Wa)=yuN,; »=0,L...,n—-1 


zur Lösung vor. In einem konvexen Bereich ® des 


BUY:::% (r—1)_ Raumes seien die exakte Lösung 
und alle betrachteten Näherungslösungen enthalten, 
und die Funktion f genüge dort für alle Punktpaare 


YYr-..4y"» und y*,y*#,...,y*@®D der Lie- 
sScHITzbedingung 


Me, 1,8. - 2,9" 9) — fe yr.9*,. u" 9)| 


n—1 
Sk 2 W"—yv*"] ee 


Die Berechnung von diskreten Näherungswerten 
der Lösung mit Differenzenschemaverfahren besteht 
aus zwei Schritten, der Berechnung des ‚‚Anfangs- 
stückes“ und der ‚‚fortlaufenden Rechnung“. Es 
werde vorausgesetzt, daß das Anfangsstück berech- 


net ist und damit Näherungswerte vr u aD Per 
nv =l,l,...,n—1 an den Stelln z,=z, 
+ uh bekannt sind. Die exakten Werte werden im 
folgenden mit 


y%z,); v=0,L,..,.n—l, 


au Yzu); REED y" —Dz,)) = F(x,) 


bezeichnet, während die zugehörigen Näherungswerte 
als 3 und f„ geschrieben werden. Die Fehlerbeträge 


der yw seien je) , also 


vr — ya = 1] EEE 


Diese Fehler hängen wesentlich von der gewählten 
Schrittweite h ab: 


Be 0) 
a (Dee a ZELLEN. 


Für die folgenden Betrachtungen sei nun stets vor- 
ausgesetzt: 


| = Max {ui —yen]} 


— eh) > 0 für h — 0 Er A 


ee: VOLL. 
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3. Das Adamssche Interpolationsverfahren 


Ist das Anfangsstück bis zur Stelle x — x, — %+rh 
berechnet, so werden weitere Näherungen ve nach 
dem Apamsschen Interpolationsverfahren ermittelt: 


n—m—1 


h? 
ae = 2 er vr ee 
0o=0 0% 


p 
ee EN Me 


mitp<sr. 


Für die exakte Lösung der Differentialgleichung (1) 
gilt: 


n—m—1 2 
re, + ı) z z RE Re) 
p 
a a nm, 9,p Fz, +1 =.) SE Be, +1 Mm. 


Der Quadraturfehler R#_ kann in der bekann- 


m,p +1 
ten Weise abgeschätzt werden durch 
p+n—m+11pz 
er ( 


X Max |FP+N | St? A 


wobei vorausgesetzt wird, daß Fe+ (x) im gesamten 
Integrationsbereich beschränkt ist. 

Subtrahiert man die Gleichungen (7) von den 
entsprechenden Gleichungen (6), setzt u — 9 z,) 
— EN), =r—9,r—p+]1,...,r,r+ 1,50 erhält 


man mit Berücksichtigung der Lirscutrzbedingung 
(2) die Fehlerungleichungen: 


n—m—1 


[7} 
lem |< Du RT |.m+ 0) 
BT e=0 g!' 
p n—1 
E —n | (») 
+ M ER ER Br > a 
o=0 v=(0 
a ed ee SE I, 


oder 
n—1 - 

| ‚| —h"FK eg i > |] 

n—m—1 

< 19 


h? | 
0 


o=0 


Ss 
IM | 
N 
in 
} —_ 
ji 
| 
a — 


p 

t oe, ü | 

= Re TR & |Pi_m, 0,p 
o= 


\ 
Fr 


u u RE a ET (10). 
Mit (5) wird daraus: 
n—1 
(m n—m gr | I | (+) 
le |—% K Ben 2. &,11 
| ur 
n—m—1 
Fe De (m--0o) 
= 52 en e, = 
o=0 DO: 
p n—i 
ö ee r | >’ ,) 
FATTE = Bel 2 &% 
a „=0 
f p+2 „(m) : 
+ h Ede . (11) 
Dieses Ungleichungssystem (11) hat die Form 
) <R Ing > 
Ar 1<Be+t3,;, 12) 
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ed, “ 
Y=-6—hlI - . Loss <A füh<JH, 


bs1 dag ER dan 


B=E+h ‚o<sbir<Bfürk<H 
a. ba Re! ban, 
| ed) ei”) si) 
2 Eu ie S et W 
an) Bas) „nm 1) 


Die Matrix X ist monoton ([2] S. 126/127); es gilt 
daher sicher 


A,yı=dertd,yı (I) 


mit +1 = &, 41 Ei 
Nun ist 
Ar, = + B-Herrt3,rı- 
Die Matrix 8 — W hat nur nichtnegative Elemente, 


daher ist 


We >UVe, für h >20 


r+1 
und somit 2% 

EN (14) 
wegen der Monotonie von \. 


Analog (13) gelten dann, wie man leicht erkennt, 
allgemein die Gleichungen 


RT MA u age Fine (15), 


BE a ED: ,+u41 S Erarı. 


Betrachten wir jetzt (15) an der festen Stellex = M 
mit M=!h,l>r; 
Ma=Bde_, +3. 


Mit AT!8 =6 gilt, wie man durch vollständige 
Induktion leicht bestätigt: 
I—r—1 
gt. + 3 an! 3%» Kl8), 
u=0 
Die Matrix 6 — X! 8 hat die Form 


G11C12 «++» Cın 


O1 gs Tee Cam 


G=C+Al. . IJmit0 <are<(Q (17), 


Cnı Cna- +» «Cnn 


da die Reihe 
UI=-CE+Ä+MUL... 


mit AU =&— W konvergiert und die Matrizen ® 
und W nur nichtnegative Elemente haben. (Die 
Schrittweiteh wird bei dem Verfahren so bemessen, 


daß alle Eigenwerte von h X dem Betrag nach <1 
sind [2]). 


Bezeichnen wir das größte Element ‚der Matrix 
&,(u=1,2,...), mit Max(C), so gilt 


Max(6)<1-+%0, 
Max (©) <(l+he)(I+nhC)<S(l+nR0)%, 


Max (6”)< (FARO) +nKOFTIS(+nKOM. 
Für alle « < 1 gilt daher 


Max (C*) < (1 + nh 0)! 
M 


— (l+nh0)" <exp (nem) ” . (8). 


'Die Vektoren A! I _u haben die Form 


0) 
u 
n ( 
Eu MEZ an Mi, u 7 > 
—ı 
a 


Mit der Bezeichnung 
Max(0/” ,) - m), „=0,,... I—r—1 
gilt dann nach (16) endlich die Abschätzung: 


el") < exp(nC M) X el) + ar +1 u .s@)) (20), 
E; m=0(,l,...n—]1 
Mit der Voraussetzung (5) folgt daraus 
dm) - (WM) >O für h>0... (el). 


Wie man nachträglich erkennt, ist die Konvergenz 
gleichmäßig für alle betrachteten I. 


4. DasAdamscheExtrapolationsverfahren 


Die Betrachtungen verlaufen ganz ähnlich, es ist 


hier A= 7! = €. Da die Matrix E von monotoner 
Art ist, so gilt analog (16) eine Gleichung der Form 


I—r—1 
a=8 te + N BB“ 


DR 
u=0 “ 


Hieraus folgt die Abschätzung 


n—1 


ei) <exp(nb M) (e”) +AaP+tıy ag > 
r=0 
m=(0,,..,n—1, 


woraus wieder die behauptete Konvergenz folgt. 
Für n = 1 erhält man (mit anderen Konstanten) die 
Abschätzung von Gaier [1]. 

Abschließend sei bemerkt, daß sich für einige 
weitere Verfahren, wie z. B. das trotz seiner nume- 
rischen Instabilität viel benutzte Verfahren der 
zentralen Differenzen, das Verfahren von E. Lix- 
DELÖF und das STÖRMERsche Extrapolationsverfahren 


die Konvergenz auf demselben Wege nachweisen 
läßt. 
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erklärt < 
ms pten u a an den ent- 
2» Fr 1 Stellen von den gegebenen 
rten. im Sinne der kleinsten Fehler- 
wenigsten abweicht [1], [2]. Dieses Ver- 
at sich auf die Randpunkte des Intervalls 
ıpt nicht, auf die den Randpunkten benach- 
enden nur unter der ' Beschränkung p= lete. 


s 


-—. Ein Be Verfahren könnte man natürlich auf 
z periodische Funktionen, bei denen keinerlei Rand- 
punkte ausgezeichnet sind, grundsätzlich ohne Be- 
E chränkungen für p anwenden. Dagegen ließe sich 
‚jedoch einwenden, daß es der Natur vieler Probleme, 
die auf periodische Funktionen führen, angemessener 
ist, einen „glatten Verlauf im Großen“ zu erstreben, 
als die gegebenen Funktionswerte in Abhängigkeit 
gewisser Nachbarpunkte zu korrigieren. So können 
Be z. B. bei trigonometrischer Interpolation der durch 
_ Ablesung aus einer graphischen Darstellung gewon- 
E nenen Funktionswerte gewisse höhere FOURIER- 
: E= Koeffizienten auftreten, die im wesentlichen die 
Schwankungen der Zeichen- bzw. der Ablesegenauig- 
keit widerspiegeln, dem betreffenden Problem jedoch 
fernliegen. Die Unterdrückung der ‚‚oberen Hälfte“ 
' aller Fourıer-Glieder, die bei der trigonometrischen 
Interpolation auftreten können, bezeichnen wir daher 
als „trigonometrische Glättung“. Diesen zunächst 
sehr grob beschriebenen Prozeß werden wir anschlie- 
ßend auf das Prinzip der kleinsten Fehlerquadrate 
zurückführen. Dies geschieht unter weitgehender 
Analogie zur eingangs beschriebenen bekannten Glät- 
tungsmethode, indem an die Stelle von Polynomen 
nunmehr FOuRIER-Polynome treten. Schließlich er- 
gibt sich eine für praktische Zwecke leicht zu hand- 

-  — habende numerische Glättungsformel. 


Eine mit 27 Der usche Funktion sei an 2 N äqui- 


 distanten Stellen 
7 
m=T BEN lh. B), 


numerisch durch f„, 2 =0...(2 N —1), gegeben 

(N =1). Diese Funktionswerte lassen sich durch 

genau ein Fourıer-Polynom mit 2 N Freiheitsgraden 
von der Form 


2 fg) = En eos up + Zu sinup ..(2) 


igonomirch interpolieren derart, daß f(p.) = fu, 

=0...(2 N—]), wird [1],[2]. Die durch tri- 

Be seksieche Glättung aus den f„ hervorgehenden 

- geglätteten Funktionswerte /„ erklären wir durch 

die (an entsprechenden Stellen genommenen) Werte 

eines Foukies-Polynoms von der halben Ordnung 

wie (2), welches an den Stellen p, von den gegebenen 

f„im Sinne der kleinsten Fehlerquadrate am wenigsten 

abweicht. Dazu erweist sich eine Fallunterscheidung 
nach geraden und ungeraden N als notwendig. 

Der Fall ungerader N ist insofern einfacher, 
weil sich ohne weiteres sagen läßt, was man unter 
einem FoukIEr-Polynom von der halben Ordnung 
wie (2) zu verstehen hat, nämlich 


A Bel, 


5 


r 


2 
Z a„cosup+ & 
u=0 1 


u=1 


ip) = businup . (8), 


a I a Aida Bat I Eu u a En un 


Sc a Dee De 


führt, dann i in bekannter a auf di die a ‚BEssEL 
i Bo Formeln“ bezeichneten Elehinape je ‚2 es: Ben 


er N Fe 
rn I N-1 
= COS U Pr ; Mena A 

“ N eh j 2 

b,= NY 2 Mayor Nino 5 
Durch Einsetzen in (3) folgt 
7 a 

1 22-1 


IQ=;y 2 te cos (pr — nd 


auf Grund der leicht zu beweisenden Identität 


Zn+1l)y 
re 
142 oup- zen _ en 
nu=—n sin Ki (4), 
2 
De) 
wird 
| 2x sin a 
ig) = fi 
2N = sin il 


Die gesuchten geglätteten Werte Ti erhält man aus 
dieser Gleichung, indem man 9> gu u=0... 
(2 N —]), streben läßt: 


Rn in PZWR 
ee 
„=glurtoy ge En burn" 

„tu Nora 


Setzt man f„ für «>2N durch f„= fu-en fort, 
so folgt 


Be 2Nn—ı 
Deohart 2 Gruss, K=0: (ZN —1)(5), 
ne 
»+gerade 
mit 
v—1 
Belt 

Bell) >, m 
2N sin ® - 

2 

y—=1...(2 N —]I), v ungerade 


Im Falle gerader N kommen die FOURIER-Poly- 
nome 
7 z 


fo) = E,arcosup+ 2 2 bisinup. (7), 
Die 
z ei 

f**@) = Kaftcosup + Zip sinup . (8), 
u=0 


mit N—1 bzw. N +1 Freiheitsgraden der Vor- 
stellung eines FouRIER-Polynoms von der halben 
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Ordnung wie (2) am nächsten. Der Fall mit N Frei- 
heitsgraden muß bemerkenswerterweise ausscheiden, 
da sonst der Ansatz nicht für alle vorkommenden 4 
gleichzeitig Glieder mit cos up und sin up enthalten 
könnte; eben dies müssen wir aber verlangen, um beim 
Glättungsverfahren keinerlei Funktionswerte vor an- 
deren auszuzeichnen. Es zeigt sich nun, daß die An- 
sätze (7) und (8) auf Glättungsformeln führen, die mit 
(5) nicht den Vorteil, daß nur Glieder für ungerade v 
auftreten, gemeinsam haben. Dies leistet dagegen das 
Mittel 


-- 1 
fo) = 2 He + FO) - 
das für unsere Zwecke offenbar völlig gleichwertig 


mit (7) bzw. (8) ist. Auf Grund der Besserschen 
Formeln 


2-1 N 
a=a*t=m 2 Jreosup R=1..(g—1 » 
v=0 
1 2N—1 { N 
y-bit-r 2 fsinup, u=|l vs ) 
ur ee n 
“ NV 2 f» cos 2 
. 1 2N—1 2» 
u = 
® N 2 f, sin 2 


vo| 


—ı 
3X | +2 ER cos u (Pp, — P) cos ven 


nu=1 
Benutzung von (4) führt auf 


” ı 29—1 
Id = 2, f, sin 8 —. 
Fe 


- 


N (9, —p) 
5 ct 


Grenzübergang 9> 9, u =0...(2? N—1), ergibt 


1 2N—1 


. W—u)n — u)n 
3N 2 ü 2 E 8 3 « 5 
7 


- 1 
fu a D In SIE 
Damit werden wir auch für gerade N auf die Glät- 


tungsformel (5) geführt, wobei allerdings die Koeffi- 
zienten @, jetzt folgendermaßen lauten: 


v—1 ) 
ee (—1) x Pr L 


ON Rn U. (10). 
v=1...(2 N—1), » ungerade 


Die verschiedene Gestalt der „‚Glättungskoeffizien- 
ten“ (6) und (10) für ungerade bzw. gerade N läßt 
sich auf folgende gemeinsame Form bringen: 


v—] 
—ı) ? R 
G, => to a BEE N I 9 
4N (etz ee) (),\ (11). 
v=1...2N—|]), » ungerade 


Dieser Darstellung entnimmt man übrigens leicht die 
Symmetrieeigenschaft 
Gy, =%, v=1... 


(2N—1), » ungerade (12), 


die von vornherein auch erwartet werden mußte. 
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Beispiel 1. Essei N > 3, und die gegebenen Funk- 


tionswerte lauten 
fu=cop, #=0...@N—). 


Die trigonometrische Interpolation (2) liefert dann 
aus Eindeutigkeitsgründen 


f(p) = eos P. 


Da es für N > 3 stets eine durch (3) bzw. (7) bis (9) 

egebene Funktion f(p) gibt, die an den Stellen 94, 
a = 0...(2 N — 1), mit f(p) übereinstimmt, nämlich 

Ip)=co9, 

so ist diese Funktion wiederum aus Eindeutigkeits- 
gründen die gesuchte Approximation im Sinne der 
kleinsten Fehlerquadrate. Es folgt daher f„ = fü 
u=0...(2 N—]1), d.h. die Durchführung der 
trigonometrischen Glättung bleibt wirkungslos. Die- 
ses Beispiel zeigt deutlich den Unterschied zum Glät- 
tungsverfahren mit Hilfe von Polynomen, welches 
die unkBonweie in diesem Falle verändern würde. 


Beispiel 2. Für N = 6 lauten die Glättungskoeffi- 
zienten gemäß (10) bzw. (11) ? 


G, = Gun — 0,3110042 . 
G,= 6, = — 0,0833333 , 
G,=G = 0,0223291. 


Für die gegebenen und die zugehörigen geglätteten 
Funktionswerte sei dann das folgende Zahlenbeispiel 
gewählt: 


AR Fu u 
0 0,7900 0,7351 
1 0,4400 0,4004 
2 0,1300 0,2086 
3 0,3000 0,2451 
4 0,3700 0,3976 
5 0,5400 0,4930 
6 0,3900 0,4832 
7 0,5800 0,4762 
b) 0,5500 0,5947 
9 0,7700 0,8265 

10 1,1500 1,0157 

1l 0,8600 0,9937 


Zur Veranschaulichung der Wirkungsweise der tri- 


gonometrischen Glättung sind die /„ und f„ im Bild 
aufgetragen und jeweils nach (2) trigonometrisch 
interpoliert. 


LEERE er 


nn 


Fa u 
) 5 X 14 er 
Trigonometrische Interpolation der gegebenen (—*k—) und 
der geglütteten (—o—) Funktionswerte 
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_ Die meisten Iterationsverfahren zur Auflösung 


I; 


we er. BES Di 
Kleine Mitteilungen 


der Iterationsverfahren 
ichungssysteme im 


eines linearen Gleichungssystems von n Gleichungen 


_ mit n Unbekannten, in Matrizenschreibweise 


u a Dre EEE CHE, (1) 


fl 
u 
5 
P 
BD 


(Gesamtschrittverfahren, Einzelschrittverfahren, 
Herrwıssches Verfahren in Gesamt- und Einzel- 
schritten, Verfahren von v. Mıszs usw.) beruhen auf 
einer Zerlegung der Matrix Y in zwei Summanden 
B und € und der Iterationsvorschrift 


SD EEE Ye ne (2). 


Dabei bedeutet +” die v-te Näherung für die Lösung x. 


_ Unter der Voraussetzung det ® + 0 folgt 
N he a a 


wobei x‘ ein beliebig wählbarer Ausgangsvektor 
sein soll. j 

Im Falle det Y = 0 ist die Lösung eindeutig be- 
stimmt und daher für die Konvergenz des Verfahrens 
notwendig und hinreichend, daß 


Im — 810 =), 

„>00 
d. h. die Eigenwerte A, von &® = — 8”! @ im Innern 
des Einheitskreises liegen. 


Ist jedoch det Y = 0, so liegt mindestens ein Eigen- 
wert auf dem Rande des Einheitskreises, denn aus 
det A=0 folgt det ($G—&)=0, d.h. — BE 
hat u. a. den Eigenwert 1. Das bedeutet jedoch nicht, 
daß obiges Verfahren versagen muß, wenn das System 
(1) überhaupt eine Lösung besitzt. Eine solche Lö- 
sung ist ja nicht eindeutig bestimmt, vielmehr wird 
der bei konvergenter Iteration (was die Existenz 


von lim (— 816)’ = fi bedingt) erhaltene Vektor _ 


,r>o0 
rt = lim x® von der Wahl des Ausgangsvektors 2” 


‚>00 
abhängen, und die Forderung 


EEE} —EN]-% 


zieht nieht die Forderung = D nach sich. 


Es gilt folgender Satz: 

Vorgelegt sei ein lineares Gleichungssystem (1) mit 
det X = 0. Das System sei lösbar. Zur Ermittlung 
der Lösung werde ein Iterationsverfahren der Form 
(2) mit det 8 = 0 benutzt. Die Matrix & = — B1C 
besitze außer den Eigenwerten, die den Wert 1 haben, 
nur solche, die dem Betrage nach kleiner sind als 1. 
Die zu den Eigenwerten 1 gehörenden Elementar- 
teiler seien sämtlich linear. Dann konvergiert das 
Verfahren und liefert eine Lösung des Systems. 


Beweis: 
Aus der Iterationsvorschrift (2), d.h. 
r” = & ıG =» u Dal Dr 
folgt 
= Er + E+G+G?+--- + GI) B-1D. 


Fall 1) v= o (homogenes System). 


Unter den Voraussetzungen des Satzes konvergiert 
die Folge {®”} gegen eine Matrix $; denn ist 9 die 
JorDAnsche Normalform von ®, also G=-TYIT 
mit einer geeigneten Transformationsmatrix T, so 
ist &= TFT, wobei % = lim %”, und diese Grenz- 


vo 
matrix existiert. Sie enthält in der Hauptdiagonale 


Een 


nur Nullen und Einsen und sonst Nullen. Somit 
liefert die Iteration einen Vektor 


r = lim & 0 = 870, 


v>o 


Ar= 0 genügt. r wird i.a. vom Nullvektor ver- 
schieden sein, wenn man r‘® + o wählt. 


Fall 2)v +. 
Da nach Voraussetzung eine Lösung x, des Sy- 


stems (1) existiert, folgt 


Brr=nHn+Bicyh= (En 


‚und daher | 
WOW HE-Nu-® W+n- En 


Nach Fall 1) streben hier das erste und dritte Glied 
der rechten Seite jeweils gegen eine Lösung des homo- 
genen Systems und damit die Folge {r®} gegen eine 
Lösung des inhomogenen Systems. q.e.d. 


Natürlich ist mit dem so bestimmten EIN N 


+4 8% ‚die Forderung Kl) x] =» er- . 


füllt, denn sie führt auf die Bedingung 2? = $, was 
wiederum mit %=% erfüllt ist. Daß jedoch $% 
idempotent ist, ist trivial. 


Selbstverständlich ist es schwierig, wenn nicht 
sogar praktisch oft unmöglich, bei einem vorgelegten 
System und einem bestimmten Iterationsverfahren 
festzustellen, ob die Voraussetzungen des Satzes er- 
füllt sind oder nicht. Doch gibt es Fälle, wo diese 
Frage entschieden werden kann. 


Z. B. sind sie gewiß dann erfüllt, wenn & symme- 


trisch ist und für die Eigenwerte A, von © die Be- 
ziehung 0 <A,=1 besteht. Dieser Fall liegt vor, 


wenn Y symmetrisch und positiv semidefinit ist (was 
bekanntlich durch linksseitige Multiplikation des Sy- 


‚stems (1) mit J{’ erreicht werden kann), und man zur 


Lösung des Systems (1) das Verfahren von v. Mıszs [1] 
benutzt. . Bei diesem ist & = & + d mit geeignetem 
d. Auf die Tatsache, daß speziell dieses Vorgehen 
stets zu einer Lösung führt, wies bereitsQUADE [2]hin. 


Die Voraussetzungen des Satzes sind ferner bei 


“ Anwendung des Einzelschrittverfahrens erfüllt, wenn 


9{ eine symmetrische, positiv semidefinite Matrix 
vom Range n — 1 ist, deren Hauptdiagonalelemente 
alle von Null verschieden (also positiv) sind. Unter 
diesen Bedingungen ist A, = 1 nur einfache Wurzel 


der charakteristischen Gleichung von & (also der 
zugehörige Elementarteiler von der Ordnung 1), 
während alle restlichen Eigenwerte von & dem Be- 
trage nach kleiner als Eins sind. Der Beweis dieser 
Tatsache kann z.B. sofort dem entsprechenden 
Beweis von ÜOLLATZ für positiv definite symmetrische 
Matrizen [3] nachgebildet werden. Ob die Forde- 
rung Re =n—1 erfüllt ist, kann häufig schon 
von der Problemstellung her entschieden werden. 
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klassischen Funktionentheorie, die ja auf einem 


W. Lietzmann, Experimentelle Geometrie. 
111 S. m. 132 Abb. Stuttgart 1959. B. G. Teubner 
Verlagsges. Preis geb. 12,60 DM. 


Bei der Durchführung von geometrischen Kon- 
struktionen mit realen Zeichengeräten hat man es 
nicht mit Elementen der Präzisionsmathematik, wie 
Punkte, Geraden und Ebenen, zu tun, die durch ein 
Axiomensystem festgelegt sind, sondern mit Dingen 
der realen Umwelt. Man stößt daher, abgesehen von 
den Mikro- und Makrokonstruktionen, die sich prak- 
tisch überhaupt nicht durchführen lassen, insofern 
auf Schwierigkeiten, da bereits jeder Punkt einer 
realen Zeichnung als räumliches Gebilde aufzufassen 
ist. So entstehen gegenüber der Präzisionsmathe- 
matik bei jeder realen Konstruktion Ungenauig- 
keiten und Fehler. Einen Ansatz zur Berechnung 
solcher Fehler wird im ersten Abschnitt vorliegenden 
Büchleins geboten. 

Sodann wird auf Konstruktionen hingewiesen, die 
sich 
a) nur mit dem Lineal (Steiwersche Konstruktionen), 


b) nur mit dem Zirkel (Mascmaronısche Konstruk- 
tionen), 

e) mit Zirkel und Lineal, 

d) mit anderen Zeichengeräten, wie z. B. mit Stech- 
zirkel, Trisektionszirkel, usw 


durchführen lassen. Dieser Teil erinnert an das 
bekannte Büchlein von HJELMSLEVY ‚‚Geometrische 
Experimente“. 

Anschließend bringt der Autor eine Fülle von 
interessanten Beispielen von Falt- und Spiegel- 
konstruktionen, bespricht sodann Scheren- und 
Fadengeometrie und Messungen überhaupt. Bei- 
spiele aus der dynamischen Geometrie, die teilweise 
auch dem Gebiete des Sports entnommen sind, be- 
schließen dieses lesenswerte Büchlein, das viel ein- 
schlägiges Unterrichtsmaterial in sehr geschickter 
Aufmachung und mit vielen Textabbildungen bietet. 


Dresden R. BEREIS 


I. N. Vekua, Systeme von Differentialglei- 
chungen erster Ordnung vom elliptischen 
Typus und Randwertaufgaben. (Mathemati- 
sche Forschungsberichte, Band II). 107 S. Berlin 1956. 
VEB Deutscher Verlag der Wissenschaften. Preis 
brosch. 10,40 DM. 


Das vorliegende schöne Buch beschäftigt sich mit 
Systemen partieller Differentialgleichungen der Form 


ou 0) 
2 | 2 
er. 2 au+rbv-+f, 
ou | vo. i NT 
&y 2 cuw+dv-+g 


(a,b,c,d,f,y gegebene bekannte reelle Funktionen), 
oder, was dasselbe bedeutet, mit der komplexen 
Differentialgleichung 


oU 


02 


=-AU+BUH+HF.....0) 


U=u+iwz=x-+iy;4,B,F bekannte kom- 
plexe Funktionen). Nach einer Einführung ($1) 
wird die Klasse der Funktionen U=u-iv mit 
stetigem dU’/dz betrachtet, und zwar unter Verwendung 
einer neuen (allgemeineren) Definition für die Ope- 
ration d/dz ($ 2). 

Die Theorie der Funktionen dieser Klasse, die 
gleichzeitig den Gleichungen (1) bzw. (2) mit stetigen 
Koeffizienten genügen (die ‚regulären Lösungen“), 
hat, wie sich dann zeigt, viel Gemeinsames mit der 
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Spezialfall von (1), den CavcHy-Rızmansschen 
Differentialgleichungen, aufbaut. So ergeben sich 
u. a. Sätze über den Zusammenhang mit Differential- 
gleichungen zweiter Ordnung (Verallgemeinerung des 
Larraceschen Operators) und über eine Verall- 
gemeinerung der CaucHYschen Integralformel, die 
sich als Lösung einer gewissen FrepHoLmschen Inte- 
gralgleichung ergibt, sowie über die gleichmäßige 
Approximation und Reihenentwicklung. _regulärer 
Lösungen ($$ 3—7). ? 
Einen breiten Raum nehmen die Randwertauf- 
gaben zu (1) bzw. (2) mit Randbedingungen der 
Gestalt zu+ ßv=y (a, f,y HöLDER-stetige Rand- 
funktionen, «+ ß? +0) bei glatten Randkurven 
ein, die auf» singuläre Integralgleichungen vom 
Normaltypus (vgl. MuscHELiscHVILL, Singuläre Inte- 
gralgleichungen) zurückgeführt und in dieser Form 
ausführlich diskutiert werden. Hier ordnet sich z. B. 
das Problem der ‚‚derivee oblique‘ bei elliptischen 
Differentialgleichungen zweiter Ordnung ein. Für 
einige näher betrachtete Spezialfälle ergeben sich 
schließlich einfachere Lösungswege ($8). Als wich- 
tiges Anwendungsgebiet wählt der Verf. die Theorie 
der elastischen Schalen. Er stellt zunächst einige 
Grundtatsachen aus der Flächen- und Schalentheorie 
zusammen und formuliert die Randwertprobleme des 
momentefreien Spannungszustandes einer Schale, die 
auf Grund der vorangehenden theoretischen Ergeb- 
nisse diskutiert werden ($9). Im Hinblick auf die 
schwierige praktische Herstellung der Lösungen von 
(1) bzw. (2) wird zum Abschluß der Fall analytischer 
Koeffizienten gesondert untersucht und gezeigt, daß 
sich die Integration der Gleichungen hier einfacher, 
nämlich durch Zurückführung auf eine VOLTERRAsche 
Integralgleichung, durchführen läßt ($ 10). 


Es ist sehr zu begrüßen, daß der vorliegende, 1952 
erschienene Forschungsbericht des gerade um dieses 
Gebiet sehr verdienten Autors jetzt auch in einer 
deutschen vorzüglichen Übersetzung vorliegt, die 
man seinem Schüler W. Schmpr verdankt. Nicht 
zuletzt muß auch das reichhaltige Literaturverzeichnis 
hervorgehoben werden, in das auch noch eine Reihe 
später erschienener Arbeiten aus diesem Gedanken- 
kreis aufgenommen wurde. 


Gießen K. MARUHN 


K. R. Biermann, Johann Peter Gustav Lejeune 
Dirichlet, Dokumentefürsein Leben und Wir- 
ken (zum 100. Todestag). (Abhandlungen der Deut- 
schen Akademie der Wissenschaften zu Berlin.) 868. 
m. 4Abb. Berlin 1959. Akademie-Verlag. Preis 
brosch. 10,50 DM. 


Die vorliegende Schrift enthält interessante Bei- 
träge zur Geschichte des Lebens und Wirkens von 
DIRICHLET, die auf Studien von Archivalien beruhen. 
Nach einem Überblick über das auf DIRICHLET bezüg- 
liche Schrifttum, wobei die Namen KUMmMER, Bor- 
CHARDT, KRONECKER, MINKOWSKI, KLEIN, LAMPE, 
VARNHAGEN und vor allem A. v. HUMBOLDT zu nennen 
sind, wird an Hand von veröffentlichten Akten zu- 
nächst über die Vorgänge bei seinem Eintritt in den 
preußischen Staatsdienst berichtet. Die Habilitation 
in Breslau, die nicht ganz reibungslos verlief, wird in 
interessanter Weise geklärt, ebenso die bald danach 
erfolgte Beurlaubung und Aufnahme der Unterrichts- 
tätigkeit in Berlin sowie seine dortige Ernennung zum 
ao-Professor. Bekanntlich war auch die Vorlesungs- 
tätigkeit in Berlin anfänglich mit Schwierigkeiten ver- 
bunden, worüber ein vollständiges Verzeichnis seiner 
Vorlesungen von 1829—55 einschließlich der Hörer- 
zahlen wertvollsten Aufschluß gibt; man weiß, welche 


du 


De 
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Bedeutung diese Vorlesungen für die Entwicklung 
des mathematischen Hochschulunterrichts in Deutsch- 
land gewonnen haben. Weitere Akten beziehen sich 
auf sein Wirken an der Allgemeinen Kriegsschule in 
Berlin, auf den bedeutenden Einfluß, den CRELLE und 


. A.v. HUMBOLDT zur gerechten Würdigung der genialen 


Persönlichkeit DIRICHLETS ausgeübt haben, auf sein 
Wirken in der Berliner Akademie, auf seine (von ihm 
abgelehnte) Berufung nach Heidelberg und auf seinen 
Weggang von Berlin nach Göttingen. Ein Anhang 
mit einer genauen Zeittafel seines Lebens und Wirkens, 
Akten aus dem Zentralarchiv (Abt. Meseburg), Ver- 
zeichnis der in Göttingen gehaltenen Vorlesungen, ein 
sehr umfangreiches Auswahl-Verzeichnis der DIRICH- 
LET-Literatur, sowie einige zeitgenössische Bilder und 
die photokopierte Wiedergabe eines Briefes an Kun- 
MER beschließen die wertvolle Schrift. 


Berlin W. SCHMEIDLER 


1.1. Priwalow, Einführung in die Funktionen- 
theorie. Teil II. (Mathematisch-Naturwiss. Biblio- 
thek, Band 22). V + 194 S. m. 28 Abb. Leipzig 1959. 
B.G. Teubner-Verlagsgesellschaft. Preis geb. 8,— DM. 


Der zweite Band dieser wertvollen Einführung 
bietet zusammen mit dem ersten im wesentlichen den 
Stoff eines einsemestrigen Anfängerkollegs über Funk- 
tionentheorie. Dabei wird an einigen Stellen, die 
dem Autor besonders am Herzen liegen, weiter aus- 
geholt und in die Tiefe gegangen. Im folgenden der 
Inhalt: 

I: Einführung des Integralbegriffes, der CAvcHY- 
sche Integralsatz und die CaucHvsche Integral- 
formel. Besonders sei erwähnt, daß hierbei auch aus- 
führlich über die Grenzwerte von Integralen vom 
CaucHyschen Typus gesprochen wird. 

II: Reihen analytischer Funktionen und Potenz- 
reihen. Entwicklung analytischer Funktionen in 
Potenzreihen. 

III: Isolierte singuläre Punkte einer eindeutigen 
Funktion. Neben der LaAurententwicklung und den 
üblichen Folgerungen daraus werden hier einige ein- 
fache hydrodynamische Anwendungen eingeschaltet. 

IV: Die Residuentheorie mit einigen Anwendungen. 

V: Der Satz vom PıcArRD, der auf Grund der Sätze 
von BrocH und LANDAT, also unter Vermeidung der 
Modulfunktionen, bewiesen wird. 

VI: Unendliche Produkte und deren Anwendung 
auf analytische. Funktionen. Hier findet sich ins- 
besondere der WEIERSTRASSsche Produktsatz sowie 
die Verallgemeinerung des Satzes über die Ein- 
deutigkeit analytischer Funktionen. Der Satz von 
MıTTAG-LEFFLER über meromorphe Funktionen wird 
nur gestreift. 

VII: Analytische Fortsetzung mit Beispielen ein- 
deutiger und mehrdeutiger Funktionen. 

Auch dieser Band eignet sich dank der vortreff- 
lichen Darstellung und Übersetzung, sowie den 
jedem Kapitel beigegebenen Übungsaufgaben sehr 
gut zum Selbststudium. 


Gießen K. MARUHN 


M. Schuler (Dr.-Ing., b. ao. Prof. em. Univ. Göt- 
tingen), Mechanische Schwingungslehre. Teil 
II: Mehrfache Schwinger. 1./2. Aufl. V + 150 S. m. 
100 Abb. Leipzig 1959. Akad. Verlagsges. Geest & 
Portig K.-G. Preis geb. 14,— DM. 


Nach einer Einführung (14 8.) in die FOURIER- 
Analyse werden in zwei Abschnitten (79 S., 47 8.) 
gekoppelte Schwingungen (Kraftkopplung, Träg- 
heitskopplung und Reibungskopplung) und Wellen- 
bewegungen (Saite, Stab und Oberflächenwellen 
idealer Flüssigkeiten) behandelt. 


Freiberg/Sa. D. RüDIGER 
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C. Berge, Espaces Topologiques — Fonctions 
Multivoques (Collection Universitaire de Math6- 
matiques) XII + 272 S. m. 47 Abb. Paris 1959. 
Dunod Editeur. Preis geb. 3.400 F. 

Das vorliegende Buch enthält manches, was man 
in einem Lehrbuch der analytischen Topologie nicht 
zu finden erwartet. So werden z.B. in einem Kapitel 
von etwa 70 Seiten eingehend konvexe Mengen und 
konvexe Funktionen im n-dimensionalen euklidischen 
Raum untersucht. Das Buch ist sehr klar geschrieben 
und setzt nur ein Minimum an Vorkenntnissen voraus. 
En Beispiele illustrieren den dargebotenen 

toff. 

Bei der großen Zahl verschiedenartigster Begriffe 
vermißt man ein Sachverzeichnis. Außerdem wird 
die Ansicht des Verlages, daß es sich hier um eine 
einfache und zugleich aber auch vollständige Dar- 
stellung des Gebietes handelt, vom Referenten in bezug 
auf die zweite Angabe nicht ganz geteilt. Eine Ein- 
führung in die allgemeine Topologie sollte z.B. auch 
die Vervollständigung eines metrischen Raumes und 
den mit den Urysounschen Metrisationssätzen in 
Verbindung stehenden Fragenkomplex behandeln. 
Wenn das Buch schon auf die fundamentalen Sätze 
über BanAacHräume eingehen will, so darf man den 
gerade für verschiedene Anwendungen sehr wichtigen 
Banachschen Satz von der inversen Abbildung nicht 
übersehen. Der auf S. 70 formulierte Satz, nach dem 
die Unterräume eines beliebigen normalen Raumes 
auch alle normal sein sollen, ist (bekanntlich) falsch. 

Wir gehen ganz kurz auf die einzelnen Kapitel ein. 
Nach drei einführenden Abschnitten (Mengen, Filter- 
basen, Abbildungen, geordnete Mengen, Zorxsches 
Lemma usw.) werden in drei weiteren Abschnitten die 
grundlegenden Sätze und Definitionen über topolo- 
gische Räume angegeben (topologische Räume, stetige 
Abbildungen, reguläre und normale Räume, kompakte 
Mengen,- zusammenhängende Mengen, topologische 
Produkte, metrische Räume u.a.). Kap. VII be- 
schäftigt sich mit Vektorräumen (lineare Mannig- 
faltigkeiten, Kegel, konvexe Mengen, Trennungssatz 
von KAXKUTANI, Satz von Haun-BAnAcH usw.), 
während der folgende Abschnitt, wie schon eingangs 
erwähnt wurde, den konvexen Mengen und Funk- 
tionen im R*R gewidmet ist. Hier findet man auch den 
Brovwesschen Fixpunktsatz und einen Paragraphen 
über bistochastische Matrizen. Das inhaltsreiche 
Buch schließt mit einem Kapitel über lineare topolo- 
gische Räume, das u. a. auch die wichtigen Fixpunkt- 
sätze von SCHAUDER und TYcHoxoFF enthält. 


Dresden M. LANDSBERG 


H. F. Schwenkhagen, Allgemeine Wechsel- 
stromlehre. Zweiter Band: Vierpole, Leitungen, 
Wellen. XII + 441 8. m. 335 Abb. Berlin/Göttingen/ 
Heidelberg 1959. Springer-Verlag. Preis geb. 39,—. 


Die „Allgemeine Wechselstromlehre“ ist 
nunmehr durch das Erscheinen des zweiten Bandes 
zum Abschluß gekommen. Der inzwischen leider ver- 
storbene Verf. verfolgte bei der Abfassung seines Wer- 
kes das Ziel, die in ähnlichen Lehrbüchern herkömm- 
liche Trennung der Wechselstromlehre in die Belange 
der Energie- bzw. Nachrichtentechnik zu vermeiden, 
ausgehend von der Feststellung, daß die Grundlagen 
der Wechselstromlehre unabhängig von der Frequenz 
sind. Er geht dabei in konsequenter Weise so vor, 
daß zunächst für irgendeinen Komplex die allgemein 
gültigen Grundlagen entwickelt und dann erst daraus 
die Spezialfälle abgeleitet werden, die dann zu den 
speziellen Anwendungen in der Starkstrom- bzw. 
Schwachstromtechnik führen. 

Der erste bereits 1951 erschienene Band behandelt 
ausschließlich quasistationäre Vorgänge, der vorlie- 
gende zweite dagegen ist den Ausbreitungserschei- 
nungen gewidmet. Im Einzelnen besteht er demgemäß 
aus den folgenden Hauptteilen: Vierpoltheorie, lange 
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Leitungen im stationären Zustand, Vierpolketten, 
Schaltvorgänge auf Leitungen, Wanderwellen und 
schließlich Wellenausbreitung im Raume. Für ein 
nutzbringendes Studium des zweiten Bandes ist eine 
gute Kenntnis der im ersten Band gebrachten theore- 
tischen Grundlagen erforderlich, auch sind die Ansprü- 
che, die in mathematischer Hinsicht an den Leser ge- 
stellt werden, merklich höher als im ersten Band, 
dürften aber von einem Studenten der Elektrotechnik 
nach Absolvierung des 5. Semesters erfüllt werden. 
Das Prinzip des Nebeneinanderstellens von Schwach- 
und Starkstrombelangen wird auch in diesem Band 
genau so konsequent durchgeführt wie im ersten. Bei- 
spielsweise werden gleich zu Beginn bei der Behand- 
lung der Zweipole Ortskurve und Ersatzschema für 
eine mit Vordrossel betriebene Leuchtröhre, eine 
Halbwellen-Dipolantenne, eine Drosselspule mit und 
ohne Eisen und einen belasteten Asynchronmotor 
nebeneinander diskutiert. Desgleichen beschränken 
sich in der Leitungstheorie die Betrachtungen nicht 
nur auf den Fall der Nachrichtentechnik, sondern 
auch die hochbelastete Drehstromleitung wird vom 
gleichen Gesichtspunkte aus behandelt. Solche Bei- 
spiele lassen sich beliebig vermehren. Sie geben der 
Lektüre des Buches einen ganz besonderen Reiz, der 
darin besteht, daß der Leser fortgesetzt auf Zusam- 
menhänge hingewiesen wird, die ihm fremd erscheinen 
mögen oder sogar gänzlich unbekannt sind, da er als 
Spezialist oftmals nicht über die nötige Zeit verfügt, 
über die Grenzen seines Spezialfaches hinwegzusehen 
und Beziehungen zu Nachbargebieten anzuknüpfen. 
Hinzu kommt die unvergleichlich prägnante Sprache 
und die geschickte Anlage und Einteilung des riesigen 
Stoffes, alles Vorzüge, die den erfahrenen Pädagogen 
verraten, der übrigens noch der Schule ORLIcH ent- 
stammt. 

Das Buch wird älteren Studenten (etwa vom 5. 
Semester an) großen Nutzen bringen, außerdem aber 
auch fertigen, bereits in der Praxis stehenden Inge- 
nieuren, die bei ihrem derzeitigen Wissensstande nicht 
haltmachen wollen, sondern das Bestreben haben, 
weiter zulernen. Der Verf. hat sich mit diesem seinem 
Abschlußwerk ein Denkmal geschaffen, das ihm auf 
lange Zeit die Hochachtung der Fachwelt erhalten 
wird. 


Dresden G. MIERDEL 


Handbuch der Physik, Band III/2. Prinzipien 
der Thermodynamik und Statistik. VII + 678 S. m. 
25 Abb. Berlin/Göttingen/Heidelberg 1959. Springer- 
Verlag. Preis geb. 160,— DM. 

Dieser Band des Handbuchs der Physik über 
die Prinzipien der Thermodynamik und Statistik 
beginnt mit einer Darstellung von E. A. GUGGENHEIM 
(118 Seiten, englisch) der klassischen Gebiete von 
phänomenologischer und statistischer Thermodyna- 
mik. Nach einleitenden Kapiteln über die Grund- 
lagen und Anwendungen der klassischen Thermo- 
dynamik wird die statistische Mechanik dargestellt 
und auf Moleküle, ideale Gase und andere Gesamt- 
heiten angewendet. Ein Kapitel über äußere Felder 
und eins über die historische Entwicklung schließen 
diesen sehr souverän geschriebenen Teil ab. 

Der nächste Teil über Axiomatik der T'hermodyna- 
mik von G. FALK und H. Jung (57 Seiten) nimmt als 
Originalbeitrag eine Sonderstellung ein. Die Verfasser 
charakterisieren ihn selbst: „Obwohl die eleganteste 
Form der axiomatischen Darstellung die sein mag, 
Herkunft und Motivierung der Axiome zu unter- 
drücken, haben wir uns bemüht, die zu den. Axiomen 
führenden analysierenden Betrachtungen in den 
wichtigsten Punkten darzulegen.‘‘ Diese Situation 
zeigt deutlich genug, wie wenig eine Axiomatik der 
Darstellung physikalischer Probleme angepaßt ist. 
Deshalb verwundert es nicht, daß selten mehrere 
Axiomatiker den gleichen Gegenstand in die gleiche 
axiomatische Form bringen. 


Einen gewissen Schwerpunkt des Buches stellt der 
Beitrag von A. Münster über die Prinzipien der stati- 
stischen Mechanik dar. Der Verfasser betrachtet 
dabei die Ableitung makroskopischer Eigenschaften 
aus der atomistischen Struktur als sein Hauptanliegen. 
Dabei steht die gedanklich strenge Grundlegung der 
Theorie und die Darstellung der allgemeinen Methoden 
im Vordergrund. Anwendungen werden nur gegeben, 
soweit zum Verständnis der allgemeinen Methoden 
erforderlich, weil sie in anderen Teilen des Handbuchs 
ihre ausführliche Darstellung finden. 

Nach den Hauptkapiteln über klassische Statistik 
undQuantenstatistik folgt ein Beitrag zur Begründung 
der Thermodynamik, der außerdem noch die Schwan- 
kungstheorie, Phasenumwandlungen und Stabilitäts- 
bedingungen beschreibt. „Die Thermodynamik der 
irreversiblen Prozesse‘ von J. Meıxner (111 Seiten) 
kann als Zusammenfassung der umfangreich vorhan- 
denen Buchliteratur dieses Gebiets angesehen werden. 
Es wird hauptsächlich die makroskopische Theorie 
der Zustände und Zustandsänderungen in der konti- 
nuierlich ausgebreiteten Materie untersucht unter der 
Voraussetzung, daß thermodynamisches Gleichge- 
wicht in mikroskopischen Bereichen herrscht. Die 
statistische Theorie irreversibler Prozesse wird vor- 
wiegend nur zur Begründung der ONSAGER-CASIMIR- 
schen Reziprozitätsbeziehungen herangezogen. Die 
übrige statistische Theorie ist noch zu sehr im Fluß, 
als daß eine Darstellung im Handbuch bereits gelohnt 
hätte. 

Der letzte Teil über Wahrscheinlichkeit und sto- 
chastische Prozesse von A. RAMAKRISHNAN (141 Sei- 
ten, englisch) liefert das Fundament für Wahrschein- 
lichkeitsüberlegungen innerhalb der Physik. Er wird 
vom Verfasser mit dem Ziel dargestellt, dem Physiker 
die Ergebnisse des Mathematikers auf diesem Gebiet 
in übersichtlicher und für ihn geeigneter Weise zu ver- 
mitteln. 


. Dresden W. MackE 


H. Sechilt, Elektrizitätslehre. 216 S. m. 187 
Abb. Basel/Stuttgart 1959. Birkhäuser-Verlag. Preis 
geb. 24,— DM. 

Das Buch ist eine Übersetzung aus dem französi- 
schen, die vom Verf. selbst bearbeitet wurde. Es 
bringt eine knappe und klare Darstellung der Theorie 
elektromagnetischer Felder, die an den Leser keine 
großen mathematischen oder gedanklichen Anforde- 
rungen stellt und daher mit Recht als elementar be- 
zeichnet werden kann. Die altgewohnte Einteilung 
in statische, stationäre, langsam und rasch veränder- 
liche Felder wird beibehalten. Neben den theoreti- 
schen Überlegungen findet sich eine Menge Hinweise 
auf praktische Anwendungen, wenn auch zuweilen 
nur durch Stichworte angedeutet, wie z. B. Zyklotron, 
Massenspektrograph. Durch intruktive kurze Rechen- 
beispiele, die jeweils an Ort und Stelle durchgerechnet 
werden, wird dem Lernenden ein tieferes Eindringen 
in den Stoff wirkungsvoll erleichtert. 

Es sind im Wesentlichen zwei Züge, die das vorlie- 
gende Buch aus der großen Anzahl der schon vorhan- 
denen und z. T. auch hervorragenden Darstellungen 
der KElektrizitätslehre herausheben: Einmal das 
Hauptprinzip, das in theoretischen Werken meist 
geübte und vielleicht schon selbstverständlich gewor- 
dene Verfahren, die Darstellung konsequent auf den 
Maxwertschen Feldgleichungen aufzubauen, nun 
auch für ein elementar angelegtes Werk zu überneh- 
men. Der zweite Punkt betrifft die unbedingte Be- 
vorzugung der einer direkten Messung zugänglichen 
Größen gegenüber den mehr formal eingeführten 
Hilfsgrößen. Es besitzen daher die elektrische Feld- 
stärke und die magnetische Induktion durchaus den 
Vorrang gegenüber der Verschiebung bzw. der magne- 
tischen Feldstärke. Deshalb wird z. B. der Gausssche 
Satz nicht, wie in der Elektrotechnik üblich, auf den 
Verschiebungsfluß angewendet, sondern auf den Feld- 


fluß. Damit muß aber auch die bekannte Komplika- 
tion in Kauf genommen werden, die bei der Berück- 


-  sichtigung von verschiedenen Dielektrika im Feld ent- 


steht und die Einführung des Begriffes der ‚‚freien“ 
Ladung zur Folge hat. Ähnliches gilt dann auch noch 
später für die magnetischen Feldgrößen. Hier wird 


_ nicht die Feldstärke H, sondern die Induktion B mit 


der Durchflutung in Verbindung gebracht. Die Wahl 
zwischen den beiden Möglichkeiten ist natürlich Ge- 
schmackssache, und F. Huxp hatin seinem bekannten 
Lehrbuch die Gleichwertigkeit beider Betrachtungs- 
weisen klar und eindeutig bewiesen; jedoch vom 
Standpunkt der gerade für ein Lehrbuch zu erstre- 
benden Einfachheit sollte man Begriffe wie ‚„makro- 
skopische“ und ‚‚mikroskopische‘‘ Ladungen bzw. 
Ströme vermeiden. 


Bei einer evtl. Neuauflage müßte die verfehlte Dar- 
stellung der Theorie starker Elektrolyte korrigiert 
und vielleicht auch das Sachregister erweitert werden, 
das recht knapp gehalten ist. Im übrigen stimmt Ref. 
mit dem Geleitwort W. Pauuis überein, daß das Buch 
„sowohl experimentellen und theoretischen Physikern 
als auch Elektrotechnikern aufs wärmste empfohlen 
werden kann“. 


Dresden G. MIERDEL 


Festschrift Richard Grammel zum 70. Ge- 
burtstag am 3. März 1959. (Ingenieur-Archiv, Bd. 
XXVIIL) VII + 3728. m. 217 Abb. Berlin/Göttin- 
gen/Heidelberg 1959. Springer-Verlag. Preis brosch. 
92,— DM. 


Mit einem Festband (als laufendem achtundzwan- 
zigsten Band des Jahrgangs 1959) ehrte das ‚‚In- 
genieur-Archiv‘ den 70. Geburtstag seines Begründers 
und seitherigen Herausgebers RICHARD GRAMMEL, 
und von einer stattlichen Reihe namhafter Vertreter 
seines Fachgebietes sind dazu Beiträge geliefert 
worden. 


Das Geleitwort schrieb GRAMMELSs langjähriger ver- 
trauter Freund C. B. Brezeno (Delft). Er würdigt 
hierin das große wissenschaftliche Gesamtwerk sowie 
die bedeutenden und weitgreifenden organisatorischen 
Leistungen des Jubilars und dessen vorbildliches di- 
daktisches Geschick als akademischer Lehrer; beson- 
deres Anliegen ist ihm dabei, die hohen menschlichen 
Qualitäten GRAMMELSs in den Vordergrund zu stellen, 
vor allem seine vaterländische aufopferungsvolle 
Arbeit in den ersten Nachkriegsjahren, die in der heu- 
tigen Blüte der aus ihren Trümmern 1946 wieder- 
erstandenen Stuttgarter Technischen Hochschule ihre 
Krönung gefunden hat. 


Leider ist es aus Platzmangel unmöglich, die Titel 
der 38 in dieser Festschrift enthaltenen — zum Teil 
recht eingehenden — Abhandlungen anzuführen, ge- 
schweige denn im einzelnen auf ihren Inhalt näher 
einzugehen. Die Thematik der Aufsätze überspannt 
alle aktuellen Teildisziplinen der Mechanik und ihrer 
Randgebiete, angefangen bei der gerade von GRAMMEL 
wohl mit besonderer Liebe gepflegten Kreiseltheorie, 
über lineare sowie nichtlineare Schwingungsvorgänge, 
das elastische und plastische Verhalten von Körpern 
unter mechanischen und thermischen Einflüssen (mit 
den hierbei anfallenden vielfältigen Sonderproblemen) 
bis zum Komplex der Strömungsphysik; auch die 
augenblicklich im Brennpunkt des allgemeinen In- 
teresses stehenden Weltraumgeschosse kommen in 
einer Untersuchung über die mechanischen Bedin- 
gungen beim Satellitenstart zur Sprache. Freunde, 
Kollegen und Schüler des hervorragenden Gelehrten 
haben ihm so in gemeinsamem Wirken eine Geburts- 
tagsgabe beschert, die darüber hinaus zu einem wert- 
vollen Besitz für die gesamte Fachwelt wird. 


Göttingen M. SCHÄFER 
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U. Olsson, Non-circular Bevel Gears. (Acta 
P. 256/1959, Mechanical Engineering Series, 5) XII 
+ 190 S. m. 31 Abb. Stockholm 1959. Acta poly- 
technica scandinavica publishing office. Preis brosch. 
Sw. Kr. 14.00. 


Das Buch behandelt zunächst die Berechnung der 
Wälzkurven für unrunde Kegelräder. Hierfür werden 
die Bedingungen des Wälzens allgemein für eine 
ganzzahlige Durchschnittsübersetzung, der Sonder- 
fall gleicher Drehzahlen beider Räder und als prak- 
tisch bedeutsamer Fall für gleiche Ausführung beider 
Räder ermittelt. Anschließend wird eine analytische 
Geometrie der sphärischen Kegelschnitte mit Polar- 


‚koordinaten und mit „Großkreis“-Koordinaten ent- 


wickelt, wobei die letzteren unmittelbar an die Be- 
wegungsformen der ebenen Getriebe anknüpfen, die, 
wie in der Getriebelehre bekannt ist, ein extremer 
Fall der sphärischen Getriebe sind. Umgekehrt 
können in der Ebene entwickelte Getriebe beliebiger 
Art, Rädergetriebe, Schaltgetriebe, Rastgetriebe usw. 
unter Beibehaltung ihrer Bewegungseigenschaften 
auf die Kugeloberfläche abgebildet werden, wobei 
die rechtwinkligen Koordinaten zu Großkreiskoordi- 
naten werden und die Getriebegeometrie unmittelbar 
in die von Oussox entwickelte ‚„Großkreis‘‘-Geometrie 
übergeht. In dem Buch werden Gleichungen und 
Berechnungen der Bogenlängen sphärischer Kegel- 
schnitte ausführlich dargestellt. Die Wälz- oder 
Teilkurven, die den Wälz- oder Teilkreisen runder 
Räder entsprechen, werden aus Kegelschnittelemen- 
ten zusammengesetzt und so berechnet. 


Das Buch von Ousson erhält seine besonders 
praktische Bedeutung dadurch, daß er eingehend 
die Herstellung unrunder Kegelräder mit den üblichen 
Evolventenscheibenfräsern und auf den normalen 
Fräsmaschinen zeigt, wobei die Berechnung der Teil- 
kopf- und Tischeinstellungen angegeben werden. 
Beispiele ausgeführter Unrundkegelräder zeigen dem 
praktischen Ingenieur die unmittelbare Verbindungs- 
möglichkeit sich schneidender Wellen bei veränder- 
licher Übersetzung, die sonst nur auf einem auf- 
wendingen Umweg über ebene Unrundräder und 
Kreiskegelräder verbunden werden. Da das Werk 
von Orsson im Anhang Berechnungsschemata und 
Tabellen enthält, wodurch die für die Gestaltung 
und Herstellung notwendigen Berechnungen leichter 
und sicher durchgeführt werden können, wird es 
sicherlich der Anwendung unrunder Kegelräder in 
der Praxis die Wege ebnen. 


Braunschweig R. A. Kraus 


N. N. Semenov, Some Problems of Chemical 
Kineties and Reactivity. Vol. 2. X + 168 S. 
m. 73 Abb. London/New York/Paris/Los Angeles 
1959. Pergamon Press. Preis geb. 35 s. 


Dieses offenbar für Chemiker und Technologen 
der chemischen Industrie bestimmte Buch befaßt 
sich mit der Theorie der Verbrennung und Explosion. 
Die schon früher (1927) vom Verf. eingeführte Unter- 
scheidung einer „thermischen“ und einer „ketten- 
artigen“‘ Explosion wird zum Ausgangspunkt ge- 
wählt und dabei insbesondere die Theorie der letzteren 
weiter ausgebaut. Es werden drei Grundtypen von 
Kettenreaktionen angegeben: die unverzweigte, die 
verzweigte und eine entartet-verzweigte Reaktion, 
auf die sich alle anderen zurückführen lassen sollen, 
Insofern sind die Untersuchungen im Rahmen der 
Reaktionskinetik von allgemeinem Interesse. — Die 
auftretenden Probleme werden mit Ausnahme eines 
Diffusionsvorganges auf gewöhnliche und Systeme 
von gewöhnlichen Differentialgleichungen zurück- 
geführt. Es sind sehr viele Beispiele einzelner chemi- 
scher Reaktionen behandelt und die Ergebnisse mit 
Messungen verglichen worden. 


Halle/Saale G. BERGER 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 
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432 ingegangene Bücher — Nachrichten“ 


EINGEGANGENEBÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten). se 


S. Karlin, Mathematical Methods and Theory 
in Games, Programming, and Economics. 
Vol.I+ I. X + 4338. und XI + 3868. London- 


"Paris 1959. Pergamon Press. Preis geb. each Vol. 


75 s. net. 

M. A. Naimark, Normed Rings. (Translated from 
the first russian edition.). XVI + 560 8. Groningen 
1959. P. Noordhoff N. V. Preis brosch. $ 12,—. 


G. Sansone — J. Gerretsen, Lectures on the 
Theory of Functions of a Öomplex Variable. 
Vol. 1. 4888. Groningen 1960. P. Noordhoff. Preis 
brosch. $ 12,—. 

T. Szab6, Repertorium und Übungsbuch der 
Technischen Mechanik. VII + 273 8. m. 
254 Abb. Berlin 1960. Springer-Verlag. Preis geb. 
DM 24,—. i 

D. Kappos, Strukturtheorie der Wahrschein- 
lichkeitsfelder und -räume. IV + 1368. Berlin 
1960. Springer-Verlag. Preis brosch. DM 21,80. 

J. Gazal6e, Les structures de commutation am 
valeurs et les caleulatrices num6riques. 788. 
m. 35 Abb. Paris 1959. Gauthier-Villars. Preis brosch. 
$ 3,04. 

A.0. Gelfond, The solution of equations in 
integers. 72 $. Groningen 1960. P. Noordhoff Ltd. 
Preis brosch. £ 3,75. 

D. Greenspan, Theory and Solution of Ordina- 
ry Differential Equations. VIII + 1488. New 
York 1960. Macmillan Company. Preis geb. $ 5,50. 


F.0. Ringleb, Mathematische Formelsamm- 
lung. (Sig. Göschen, Bd. 51/öla.) 7. Aufl. 3208. 
Berlin 1960. Walter de Gruyter & Co. Preis brosch. 
DM 5,80. 

A. Blaquiere, Mecanique non lineaire. Les 
oscillareurs a regimes quasi sinusoidaux. 
137 S. m. 45 Abb. Paris 1960. Gauthier-Villars. Preis 
brosch. $ 5,75. 

Jahnke-Emde-Lösch, Tafeln Höherer Funk- 
tionen. 6. Aufl. XII + 318S. m. 189 Abb. Stutt- 
gart 1960. B.G. Teubner Verlagsgesellschaft. Preis 
geb. DM 58,80. 

Ss. Vajda, An Introduction to Linear Pro- 
gramming and the Theory of Games. 768. 
London 1960. Methuen & Co. Ltd. Preis geb. 98 6 d.net. 


W. Blaschke / H. Reichardt, Einführung in die 
Differentialgeometrie. Zweite Aufl. VII-+ 
1738. m. 57 Abb. Berlin 1960. Springer-Verlag. 
Preis geb. DM 24,—. 
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H. Ziegler, Mechanik I. 3. Aufl. 2448. 
1960. Birkhäuser Verlag. Preis geb. Fr. 28,50. 
K.Schütte, Beweistheorie. XI + 355 8. Berlin 
1960. Springer-Verlag. Preis geb. DM 48,—. n 
6. Julia, El&ments TASEDE ae 
Paris 1959. Gauthier-Villars. is brosch. $ 7,93. 
F. Winckel, Technik der Magnetspeicher. XVI 
+ 6148. m. 325 Abb. Berlin 1960. Springer-Verlag. £ 
Preis geb. DM 72,—. a, 

L. Collatz, Differentialgleichungen für In- 
genieure. 2. Aufl. 1978. m. 115 Abb. Stuttgart 
1960. B.G. Teubner Verlagsgesellschaft. Preis geb. 
DM 21,60. DE. 

W. Gröbner, Die Lie-Reihen und ihre An- 
wendungen. VI + 1128. m. 5Abb. Berlin 1960. 
Deutscher Verlag der Wissenschaften. Preis geb. 
DM 22,—. | 

H. Wolf, Ausgleichungsrechnung nach der 
Methode der ehren Quadrate. 388. m. 
31 Abb. Hamburg 1960. Hanseatische Verlags- 
anstalt. Preis brosch. DM 6,80. 

H. L. Dryden, Th. von Karman and @. Kuerti, Ad- 
vances in Applied Mechanics. X + 2945. m. 
63 Abb. New York and London 1960. Academic 
Press. Preis geb. $ 9,00. 

M.A. Neumark, Lineare Differentialopera- 
toren. XII + 394 S.m.19 Abb. Berlin 1960. Akade- 
mie-Verlag. Preis geb. DM 44,—. 

K.P. Stanyukovich, Unsteady motion of con- . 
tinuous media. XV + 7458. London/Oxford/ 
Paris/New York 1960. Pergamon Press. Preis geb. 
£ 5.net. } 

U. Graf u. H.J. Henning, Statistische Metho- 
den beitextilen Untersuchungen. 3. berichtister 
Neudruck. XII + 291 S. m. 71 Abb. Berlin/Göttin- 
gen/Heidelberg 1960. Springer-Verlag. Preis geb. 
DM 37,50. 

S. Vianelli, Prontuari per calcoli statistici. 
XV + 15438. Palermo/Roma 1959. Abbacos.r. 1. — 
Editore. Preis geb. Lire 16.000. 

Grundzüge der Mathematik. Band II: Geo- 
metrie. Herausgegeben von H. Behnke und 
F. Bachmann. 6468. m. Abb. Göttingen 1960. 
Vandenhoeck & Ruprecht. Preis geb. DM 58,—. 

Kai Lai Chung, Markov Chains with Stationary 
Transition Probabilities. (Die Grundlehren der 
mathematischen Wissenschaften, Band 104.) X + 
278 S. Berlin/Göttingen/Heidelberg 1960. Springer- 
Verlag. Preis geb. DM 65,60. 
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NACHRICHTEN 


Herrn Prof. Dr. Dr. h. c. RıcHARD GRAMMEL hat der 
Verein Deutscher Ingenieure als seine höchste Aus- 
zeichnung die Grashof-Denkmünze verliehen. Fast 
gleichzeitig erhielt er von der American Society of 
Mechanical Engineers (ASME) die Timoshenko-Me- 
daille. H. 


Die Deutsche Gesellschaft für Versicherungsmathe- 
matik (vertreten durch Herrn Prof. Dr.-Ing. H. Krak- 
KE, Köln), die Gesellschaft für Angewandte Mathe- 
matik und Mechanik (vertreten durch Herrn Prof. Dr. 
H. GöRTLER, Freiburg) und die Wissenschaftliche Ge- 
sellschaft für Luftfahrt (vertreten durch Herrn Dr. 
W. SCHULZ, Braunschweig) haben einen Deutschen 


Ausschuß für Unternehmensforschung (DAUF) ins 
Leben gerufen, der sich folgende Aufgaben zum Ziel 
gesetzt hat: > 

l. die gegenseitige Information unter den Mitglie- 
dern zu vermitteln, 

2. Arbeiten und Veranstaltungen nach Bedarf ge- 
meinsam zu planen und durchzuführen, 

3. als umfassender nationaler Zusammenschluß die 
Vertretung der deutschen Interessen auf dem Gebiete 
der Unternehmensforschung aufinternationaler Ebene 
wahrzunehmen. 


. Weitere Organisationen der deutschen Wissenschaft 
sind zum Beitritt eingeladen worden. 
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Herausgeber und Hauptschriftleiter: Prof. Dr.-Ing. habil, H. Heinrich, Dresden A 27 


7, Friedrich-Hegel-Str. 31. Verlag: Akademie-Verlag 
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Die Zeitschrift für Angewandte Mathematik und Mechanik erscheint monatlich. Bezu 


Bestellgeld. Doppelheft DM 10,—. Abbestellungen können nur bis 
folgende Quartal noch geliefert. Veröffentlicht unter Lizenzuummer 
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| i Mathematische Elastizitätstheorie 
- der ebenen Probleme 


Übersetzung aus dem Tschechischen 


In deutscher Sprache herausgegeben von Dipl.-Math. Winfried Heinrich 


1960. XII, 478 Seiten — 207 Abb. — 1 Falttafel — gr. 8° — Ganzleinen DM 96,— 


Das Buch macht den Mathematiker mit einer wichtigen theoretisch-technischen Disziplin 


vertraut, deren Kenntnis für einige Ingenieurfächer unerläßlich ist, und zeigt andererseits 


dem Techniker, wie der Mathematiker seine Aufgaben sieht, wie er sie modifiziert, bearbei- 


tet und löst. Die De behandeln Probleme der ebenen Elastizitätstheorie, welche 


im wesentlichen Probleme der Lösung der biharmonischen Gleichung für Funktionen zweier 
Veränderlicher bei verschiedenen Randbedingungen sind. Sie stellen sich das Ziel, die 
Theorie des biharmonischen Problems und die Methoden zu dessen numerischer Lösung 
zusammenzufassen, sie zu ergänzen und so einander gegenüberzustellen, daß ihre Struktur 
ehr wird und aus diesen Vergleichen das maximal mögliche für die numerische Be- 
rechnung zu gewinnen ist. 

Hervorzuheben ist der mathematisch strenge, systematische Aufbau des Werkes, das über 
die Analysis und die Elemente der Funktionentheorie hinaus keine besonderen Vorkennt- 
nisse verlangt. Was an weitgehenden mathematischen Hilfsmitteln benötigt wird, ist in 
dem umfangreichen Anhang zu finden. Die praktische Anwendung der Methoden wird an 
verschiedenen Beispielen gezeigt. In deutscher Sprache fehlte bisher eine zusammen- 


fassende Darstellung dieses Gebietes. 


Bestellungen durch eine Buchhandlung erbeten 


1960. VI,0 Sen — 16 Ange, vn 2 u 


Das Werk gibt eine . sussmenfisne Da 
Gebiet der programngsseusrtn eshensuteiiten. Der | 
"Lehrbuchs entsprechend nach didaktischen Gesichtspunkten usg, 
ne der neuesten Ergebnisse. Ausgehend von einer E Ei u un Se 
‚ werden die durch sie gegebenen Möglichkeiten A Realiierung de der : 
Operationen entwickelt. Eine ‚Erläuterung spezieller ausgewählter I er 
Einblick in die typischen Baugruppen, ihre Arbeitsweise sowie in = . 
wicklungstendenzen gewinnen. x 


Das Buch istin erster Linie für “te Studierenden der: m 


‚und technischen Fachrichtungen an unseren Hoch- und Fachschulen ' 


IN VORBEREITUNG BEFINDET SICH 
BAND 2 


Prof. Dr. HELMUT WINKLER 


Elektronische Analogieanlagen 


etwa 296 Seiten — 172 Abbildungen — gr. 8° — Lederin etwa 36,— DM 
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